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Simulating  the  Management  of  Even-Aged  Timber  Stands 


Clifford  A,  Myers 


Introduction 

Good  decisions  that  lead  to  efficient,  low- 
cost  operations  are  as  necessary  for  timber 
management  as  for  any  other  business  activity. 
Gathering  of  information  on  which  to  base  such 
decisions  is  unusually  difficult  because  of  the 
peculiar  characteristics  of  the  business  of  grow- 
ing wood.  The  time  needed  to  mature  a  crop 
is  long,  and,  with  many  tree  species,  results  of 
changes  appear  slowly.  Decisions  must  be  made 
for  large  areas  of  land  and  for  numerous  pos- 
sible combinations  of  biological  and  economic 
conditions.  A  procedure  whereby  a  manager 
can  examine  quickly  the  probably  outcomes  of 
many  variations  in  management  should  be  of 
value.  Such  a  procedure  appears  especially 
useful  when  the  computation  time  required  is 
about  10  seconds  for  each  combination  of  pos- 
sible alternatives  examined.  This  approach  to 
managerial  decisionmaking  is  made  possible 
by  electronic  computers  and  the  technique  of 
simulation. 

Simulation  involves  the  creation  and  opera- 
tion of  a  model  that  is  similar  in  relationships 
to  the  natural  system  studied.  The  system  may 
be  a  factory,  a  forest,  or  any  one  of  countless 
other  subjects  of  interest.  Mathematical  simu- 
lation, where  the  model  is  a  set  of  equations, 
is  a  powerful  tool  of  management  and  investiga- 
tion. Mathematical  models  are  used  to  discover 
new  facts  and  to  test  alternatives.  Time  can 
be  compressed  so  a  sequence  of  events  can  be 
studied  in  a  fraction  of  the  time  required  for 
operation  of  the  natural  system,  thereby  permit- 
ting projections  into  the  future (Choraf as  1965). 

Manipulation  of  a  model  by  simulation  is 
especially  useful  where:  (1)  the  actual  system 
should  not  be  changed  until  decisions  have 
been  made;  (2)  possible  changes  are  numerous 
and  are  costly  in  money,  effort,  and  time;  and 
(3)  exact  data  cannot  be  obtained  because  of 
cost  or  complexity  of  the  relationships  involved 
(Churchman  et  al.  1957,  Dane  1966). 


Several  authors  have  presented  the  prin- 
ciples of  simulation  to  foresters,  and  have  shown 
how  forestry  simulators  may  be  constructed. 
These  authors  have  discussed  the  development 
of  models  (Chappelle  1966)  and  the  systems 
approach  to  forestry  problems  (O'Regan  et  al. 
1966).  Others  have  described  simulators  of 
timber  production  activities,  and  presented  in- 
formation obtained  from  their  operation  (Clut- 
ter and  Bamping  1966,  Gould  and  O'Regan 
1965).  Thompson's  (1966)  discussion  of  the 
need  for  tests  of  alternatives  in  forest  regulation 
points  out  a  value  of  simulation  for  research. 
Validity  of  many  long-accepted  concepts  of 
timber  management  can  be  tested  quickly, 
cheaply,  and  for  a  wide  variety  of  conditions. 

Gould  and  O'Regan  (1965)  performed  a 
special  service  to  timber  managers  by  publishing 
a  simulator  at  a  time  when  examples  were 
relatively  inaccessible  to  most  foresters.  Their 
computer  program  will  long  serve  as  a  source 
of  ideas  and  procedures;  it  has  so  served  for 
the  program  presented  here. 

Program  MANAGD  (appendix  1)  is  a  means 
for  simulation  of  the  management  of  even-aged 
stands  for  roundwood  and  saw  logs.  It  contains 
provisions  for  stand  growth,  thinnings,  harvest 
cuts,  planting  of  nonstocked  areas,  and  other 
changes  in  forest  conditions.  Inputs  to  the 
program  permit  wide  choice  in  the  management 
alternatives  and  stand  conditions  to  be  ex- 
amined. Possible  options  and  alternatives  are 
presented  in  the  appropriate  sections  of  the 
program  description. 

The  program  was  written  in  A.S.A.  Fortran 
IV  (X3.9-1966)  and  tested  on  a  CDC  6400  compu- 
ter. 2  Program  organization  will  permit  ready 
modification  to  fit  local  tree  species  and  com- 
puting equipment.  Once  the  program  is  in  use, 
improved  or  additional  data  and  procedures  can 
be  inserted  easily,  as  they  become  available. 

1  Trade  names  and  company  names  are  used 
for  the  benefit  of  the  reader s and  do  not  imply 
endorsement  or  preferential  treatment  by  the 
U.  S.  Department  of  Agriculture. 


-  1  - 


Description  of  Program  MANAGD 

Program  MANAGD  consists  of  a  main  pro- 
gram and  11  subroutine  subprograms  (appendix 
1).  Content  and  purpose  of  each  routine  are 
given  in  the  sections  that  follow.  Definitions 
of  variable  names  are  listed  with  the  source 
program  in  appendix  1  and  in  the  listing  of 
contents  of  the  data  deck.  The  test  problem 
described  on  page  14  and  reported  in  appendix 
2  provides  additional  explanation  of  the  program. 

The  terms  batch,  test,  and  game  are  used  to 
identify  individual  simulations  performed  with 
various  groupings  of  alternatives.  The  BATCH 
name  identifies  the  entire  group  of  tests  and 
games  to  be  completed  as  a  single  job  by  a 
computer.  A  test  consists  of  one  or  more 
games,  all  of  which  are  based  on  a  single  yield 
table  and  one  set  of  stumpage  prices.  Games 
of  a  test  differ  from  one  another  in  distribution 
of  acres  by  age  classes,  area  planted,  and  limi- 
tations on  the  annual  cut. 


Main  Program 

The  main  program  calls  10  subroutines  in 
proper  sequence,  and  uses  counter  UK  to  call 
an  eleventh  subroutine  (OUTPUT1)  at  specified 
intervals.  The  first  four  subroutines  called 
read  the  data  deck,  compute  and  print  a  yield 
table,  and  compute  and  print  potential  volumes 
per  acre  at  each  year  of  stand  age.  A  fifth 
subroutine  generates  a  working  circle  with  the 
specified  number  of  acres  in  each  age  class. 
Subroutines  YEARS  and  ANUAL  are  the 
dynamic  components  of  the  program.  They 
provide  for  stand  growth,  thinnings,  harvests, 
losses,  and  other  changes  in  volume  and  value. 
Four  subroutines  produce  a  record  of  results. 
Numbers  of  acres  in  each  age  class  are  printed 
at  the  end  of  the  first  year  of  a  game  and  at 
the  end  of  each  decade.  Volume  and  dollar 
values  are  printed  at  the  end  of  each  game. 
An  optional  subroutine  (SUMRY)  prints  selected 
values  from  several  games  of  a  test  together 
to  simplify  comparisons. 

The  main  program  enters  BATCH  name  and 
the  number  of  tests  in  a  batch  into  computer 
memory. 


Subroutine  INPUTl 

INPUTl  is  called  once  each  test  to  read 
values  that  apply  to  all  games  of  a  test.  Values 


entered  include  stumpage  prices,  minimum 
commercial  volumes,  and  items  used  to  compute 
a  yield  table.  Controls  on  the  program  are 
entered  as  number  of  games  in  a  test,  number 
of  years  in  a  game,  and  number  of  columns 
of  OUTPUT2  to  be  summarized  by  SUMRY. 

Values  read  by  subroutine  INPUTl  are 
printed  as  part  of  the  output  of  other  sub- 
routines. 


Subroutine  YIELD 

Subroutine  YIELD  reads  four  tables  and 
computes  and  prints  yield  tables  for  managed, 
even-aged  stands.  It  is  called  once  each  test 
to  produce  the  yield  table  that  will  apply  to 
all  games  of  the  test.  Values  in  each  yield 
table  reflect  prior  decisions  on  the  frequency 
and  intensity  (table  1)  of  intermediate  cuttings. 

Computations  performed  by  the  subroutine 
follow  procedures  described  in  detail  elsewhere 
(Myers  1966,  Myers  and  Godsey  1968).  Average 
stand  diameter  (d.b.h.)  and  number  of  trees 
per  acre  just  before  initial  thinning  are  esti- 
mated from  measurements  of  stands  of  suitable 
ages  and  densities.  Number  of  trees  and  d.b.h. 
are  used  to  compute  initial  and  subsequent 
basal  areas  per  acre.  Total  cubic  volumes  per 
acre  are  computed  from  basal  area,  d.b.h.,  and 
tree  height  (table  2)  by  means  of  a  stand 
volume  equation  (statements  for  TOTO  and 
TOTT).  Total  cubic  volumes  are  multiplied 
by  factors  to  convert  them  to  merchantable 
cubic  volumes  (table  3)  and  board  feet  (table 
4). 

Increase  in  d.b.h.  due  to  initial  thinning  is 
computed  by  the  iterative  process  of  FORTRAN 
statements  6  to  12,  inclusive.  The  equation 
for  PDBHE  provides  an  estimate  of  diameter 
after  thinning  when  diameter  before  thinning 
and  the  percentage  of  trees  retained  are  known 
(Myers  1966).  Successive  percentages  are  tested 
until  d.b.h.  after  thinning,  number  of  trees, 
and  basal  area  agree  with  the  desired  thinning 
intensity  entered  as  THIN  by  INPUTl  and  as 
defined  in  table  1.  Equations  for  DBHP,  such 
as  statement  7,  express  that  part  of  table  1 
where  basal  area  is  less  than  66.3  square  feet. 

Table  1  gives  residual  basal  area  after  thin- 
ning for  various  average  stand  diameters.  The 
values  represent  one  possible  series  of  densities 
that  could  be  used  to  guide  successive  thinnings 
in  a  stand.  The  series  is  labeled  growing  stock 
level  80  to  indicate  that  reserve  basal  area  is 
80  square  feet  per  acre  when  d.b.h.  after  cutting 
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Table  1. — Basal  areas  after  intermediate  cutting  in  relation  to  average  stand  diameter 
Growing  stock  level  80  for  Black  Hills  ponderosa  pine 
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Average 
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stand  d.b.h. 

area 
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after  cutting 

per 

after  cutting 

per 

after  cutting 

per 

after  cutting 

per 

(Inches) 

acre 

(Inches) 

acre 

(Inches) 

acre 

(Inches) 
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Sq.  ft. 

Sq.  ft. 

Sq.  ft. 

Sq.  ft. 

z .  u 

12.  3 

4.0 

35.  2 

a  n 
o .  u 

SA  A 

R  n 

7  9  S 
/Z.J 

•?  1 

Z  .  X 

13.3 

4.1 

36 . 4 

A  1 
D  .  X 

S7  A 

R  1 
o .  x 

71  9 
/  J  .  z 

9  9 
z .  z 

14 . 4 

4.2 

37 . 6 

A  9 

o .  z 

SR  S 
JO  .  J 

R  9 

o .  z 

71  R 
/  J  .  o 

2 . 3 

15. 5 

4.3 

38 . 7 

6 . 3 

A 

R  7 

O.J 

74.4 

2.4 

16.6 

4.4 

39.8 

6.4 

60.3 

8.4 

74.9 

9  S 
Z  .  J 

17 . 7 

4 . 5 

41 . 0 

A 

O.J 

A1  1 

R  S 

O.J 

7  S  A 
/  J .  1 

2 . 6 

18 . 9 

4 . 6 

42 . 2 

6.  6 

62 . 0 

8.6 

75.8 

2.  7 

20. 1 

4 . 7 

43.4 

6 . 7 

A9  Q 
o  z .  y 

R  7 
o .  / 

7A  1 

/O.J 

2 . 8 

21 . 3 

4 . 8 

44 . 6 

6. 8 

63 . 8 

8 . 8 

7A  7 
/  o .  / 

2.9 

22.5 

4.9 

45.7 

6.9 

64.7 

8.9 

77.1 

3 . 0 

23. 6 

5 . 0 

46. 7 

7 . 0 

65 . 5 

q  n 

77 . 5 

J  .  X 

24 . 8 

5 . 1 

47.7 

/  .  X 

66 . 2 

Q  1 
J  .  X 

77  Q 
/  / .  ? 

1  9 

J  .  Z 

26.0 

5 . 2 

AS  r 

7  9 
/  .  z 

A7  n 
o  /  •  u 

Q  9 
7  .  z 

7R  9 
/  o .  z 

J  •  J 

27 . 2 

j  .  j 

AQ  Q 

7  1 

A7  7 
o  /  .  / 

37  •  J 

7R  S 
/  o .  J 

J  .  4 

28 . 4 

5 . 4 

q 

7 . 4 

68 . 4 

9.4 

78 . 8 

3.5 

29.6 

5.5 

51.8 

7.5 

69.1 

9.5 

79.0 

3.6 

30.7 

5.6 

52.8 

7.6 

69.8 

9.6 

79.2 

3.7 

31.8 

5.7 

53.7 

7.7 

70.5 

9.7 

79.5 

3.8 

32.9 

5.8 

54.7 

7.8 

71.2 

9.8 

79.7 

3.9 

34.1 

5.9 

55.6 

7.9 

71.9 

9.9 

79.9 

10.0+  80.0 


Table  2. — Average  height  of  dominant  and  codominant  trees  at 
various  ages,  Black  Hills  ponderosa  pine 


Main 

Site  index 

class 

stand  age 

(Years) 

40 

1  50 

|  60 

70 

-  -  -  -  -  Feet 

10 

4.5 

4.5 

4.5 

4.5 

20 

9 

10 

12 

16 

30 

11 

16 

20 

26 

40 

17 

22 

28 

35 

50 

21 

28 

35 

43 

60 

26 

33 

41 

50 

70 

30 

38 

47 

56 

80 

34 

43 

52 

61 

90 

37 

47 

57 

66 

100 

40 

50 

60 

70 

110 

43 

53 

63 

74 

120 

45 

56 

66 

77 

130 

46 

59 

69 

80 

140 

48 

61 

71 

83 

150 

50 

63 

73 

86 

160 

51 

64 

75 

88 

170 

52 

65 

77 

90 

180 

53 

66 

78 

91 
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Table  3. — Factors  for  conversion  of  stand  volumes  in  total  cubic  feet  to 
merchantable  cubic  feet  per  acre,     Black  Hills  ponderosa  pine 


Average 

Average 

Average 

stand 

Ratio 

of 

stand 

Ratio 

of 

stand 

Ratio  of 

diameter 

merchantable 

diameter 

merchantable 

diameter 

merchantable 

(Inches) 

to  total 

volume  ! 

(Inches) 

Co  total 

volume  1 

(Inches) 

to  total  volume 

5.0 

0.332 

8.1 

0.849 

11.9 

0.940 

5.1 

.355 

8.2 

.856 

12.1 

.941 

5.2 

.377 

8.3 

.  862 

12.4 

.942 

5.3 

.400 

8.4 

.868 

12.7 

.943 

5.4 

.422 

8.5 

.872 

12.9 

.944 

5.5 

.444 

8.6 

.876 

13.  1 

.945 

5.6 

.465 

8.7 

.880 

13.3 

.946 

5.7 

.  487 

8.8 

.  884 

13.5 

.947 

5.8 

.508 

8.9 

.888 

13.7 

.948 

5.9 

.530 

9.0 

.892 

13.9 

.949 

6.0 

.552 

9.1 

.896 

14.2 

.950 

6.1 

.575 

9.2 

.899 

14.4 

.951 

6.2 

.  jy  1 

9.3 

QD? 

.  VKJi. 

14.7 

.952 

6.3 

.618 

9.4 

.906 

14.9 

.953 

6.4 

.639 

9.5 

.910 

15.2 

.954 

6.5 

.659 

9.6 

.913 

15.4 

.955 

6.6 

.678 

9.7 

.916 

15.8 

.956 

6.7 

9.8 

.  yzu 

16.3 

.957 

6.8 

.710 

9.9 

.923 

16.8 

.958 

6.9 

.725 

10.0 

.926 

17.3 

.959 

7.0 

.740 

10.1 

.928 

17.8 

.960 

7.1 

.753 

10.2 

.930 

18.3 

.961 

7.2 

.766 

10.3 

.931 

18.8 

.962 

7.3 

.778 

10.4 

.932 

19.3 

.963 

7.4 

.789 

10.5 

.933 

19.8 

.964 

7.5 

.799 

10.7 

.934 

20.3 

.965 

7.6 

.809 

10.9 

.935 

20.9 

.966 

7.7 

.818 

11.  1 

.936 

21.7 

.967 

7.8 

.826 

11.3 

.937 

22.5 

.968 

7.9 

.834 

11.5 

.938 

23.3 

.969 

8.0 

.842 

11.7 

.939 

23.9 

.969 

1  To  4.0-inch  top  in  trees  6.0  inches  d.b.h.   and  larger. 

Factor  for  an  unlisted  diameter  equals  factor  for  next  smaller  listed 
diameter.     For  example,  factor  for  15.6  inches  is  .955." 
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Table  4. — Factors  for  conversion  of  stand  volumes  in  total  cubic  feet  to  board  feet 
Scribner  rule  per  acre,1  Black  Hills  ponderosa  pine 
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Ratio  of 

Average 

Ratio  of 

1 

Average 

Ratio  of 

Average 
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board  feet 
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board  feet 

stand 

board  feet 
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to 
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to 
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12 

.2 
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16 

1 
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5.46 
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1 
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.2 
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5.47 
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1 

.13 
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1.20 
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ZU.  J 
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1 

.27 

12 

6 

3.80 

16 

5 

4 . 88 
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1 

.34 

12 

7 

3.84 

16 

6 

4.90 

20.5 

5.52 

8.9 

1 

41 

12 

8 

3.88 

16 

7 

4.92 

20.6 

5.53 

9.0 

1 

48 

12 

9 

3. 91 

16 

8 

4 . 94 

9D  7 

J.  J4 

9.1 

1 

55 

13 

0 

3.  95 

16 

9 

4.96 

20.8 

J.  JO 

9.2 

1 

62 

13 

1 

3.  98 

17 

0 

4 . 98 

20.  9 

5.57 

9.3 

1 

68 

13 

2 

4.02 

17 

1 

5.00 

21.0 

5.58 

9.4 

1 

75 

13 

3 

4.05 

17 

2 

5.02 

21.1 

5.59 

9.5 

1 

82 

13 

4 

4.08 

17 

3 

5.04 

21.  2 

5 . 60 

9.6 

1 

89 

13 

5 

4. 11 

17 

4 

5.  06 

21.  3 

5 . 61 

9.7 

1 

96 

13 

6 

4.14 

17 

5 

5.  08 

21.4 

5.  62 

9.8 

2 

03 

13 

7 

4.17 

17 

6 

5.10 

21.5 

5.63 

9.9 

2 

10 

13 

8 

4.20 

17. 

7 

5.12 

21.6 

5.64 

10.0 

2. 

17 

13 

9 

4.  23 

17. 

8 

5 . 13 

21.  7 

5 . 65 

10. 1 

2. 

24 

14. 

0 

4.  25 

17. 

9 

5 . 15 

21.8 

5 .  66 

10.2 

2. 

31 

14. 

1 

4. 28 

18. 

0 

5. 17 

21.9 

5 . 67 

10.3 

2. 

38 

14. 

2 

4.31 

18. 

1 

5.19 

22.0 

5.68 

10.4 

2. 

45 

14. 

3 

4.34 

18. 

2 

5.21 

22.1 

5.69 

10.5 

2. 

52 

14. 

4 

4.  37 

18. 

3 

5.22 

22.  2 

5.  70 

10.  6 

2. 

59 

14. 

5 

4.  39 

18. 

4 

5.24 

22.3 

5.71 

10.7 

2. 

65 

14. 

6 

4.42 

18. 

5 

5.  26 

22.4 

5 . 72 

10.8 

2. 

72 

14.7 

4.44 

18. 

6 

5.27 

22.5 

5.73 

10.9 

2. 

79 

14. 

8 

4.47 

18. 

7 

5.29 

22.6 

5.74 

11.  0 

2. 

86 

14. 

9 

4 . 49 

18. 

8 

5 . 30 

22.7 

5 . 75 

11. 1 

2. 

93 

15. 

0 

4 . 52 

18.9 

5.  32 

22.8 

5.76 

11.2 

3. 

00 

15. 

1 

4.54 

19. 

0 

5.33 

22.9 

5.77 

11.3 

3. 

07 

15.2 

4.56 

19.1 

5.35 

23.0 

5.78 

11.4 

3. 

14 

15.3 

4.58 

19.2 

5.36 

23.1 

5.79 

11.5 

3. 

21 

15. 

4 

4.61 

19. 

3 

5.37 

23.2 

5.80 

11.6 

3. 

28 

15. 

5 

4.64 

19. 

4 

5.39 

23.3 

5.81 

11.7 

3. 

35 

15. 

6 

4.66 

19. 

5 

5.40 

23.4 

5.82 

11.8 

3. 

42 

15. 

7 

4.68 

19.6 

5.41 

23.5 

5.83 

To  8- inch  top  in  trees  10.0  inches  d.b.h.  and  larger. 
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is  10  inches  or  larger.  Other  stocking  levels 
are  named  the  same  way.  For  example, 
stocking  level  100  means  that  reserve  basal 
area  will  be  100  square  feet  when  d.b.h.  after 
cutting  is  10  inches  or  larger.  Basal  areas 
for  level  100  and  for  diameters  smaller  than 
10  inches  are  obtained  by  multiplying  each 
basal  area  of  level  80  by  the  amount  100/80. 
Values  for  other  stocking  levels,  perhaps  from 
50  to  160,  are  computed  similarly. 

Periodic  increases  in  d.b.h.  due  to  tree 
growth  are  estimated  by  the  equation  for  DBHO 
in  the  loop  headed  DO  23.  The  equation  used 
in  the  example  is  for  Black  Hills  ponderosa 
pine,  and  for  a  projection  period  of  10  years. 
Equations  for  other  species  or  projection  periods 
may  be  inserted  as  desired.  Intervals  between 
intermediate  cuttings  are  one  or  more  projection 
periods  long. 

It  is  often  desirable  to  make  simulations 
more  realistic  through  introduction  of  variability 
in  values  estimated  by  equations  or  contained 
in  tabulations.  For  example,  repeated  compu- 
tations of  DBHO  without  change  in  values  of 
the  independent  variables  will  always  give  the 
same  numerical  result.  In  reality,  actual  and 
estimated  values  frequently  differ.  A  way  of 
providing  variability  in  estimates  of  DBHO  is 
contained  in  the  program  segment  between 
statements  100  and  110.  Similar  statements 
could  be  written  for  other  variables. 

Variability  is  obtained  in  three  steps;  (1) 
generation  of  a  pseudorandom  number,  (2)  use 
of  this  number  as  an  independent  variable  to 
compute  the  value  of  a  residual  (range:  -0.3 
to  +0.3  inch),  and  (3)  addition  of  the  residual 
to  the  computed  value  of  DBHO.  The  pseudo- 
random number  generator,  statement  100,  is  of 
the  form: 

Xi  e  AXi_1  +  C  (modulo  M) 

(Greenberger  1961).  Values  of  all  elements  of 
the  generator  are  specified  except  for  X  s  _  j  , 
which  is  read  in  as  variable  GNTR.  The  state- 
ment to  compute  RES  is  an  empirical  distribu- 
tion function  obtained  by  fitting  a  polynomial 
to  the  normally  distributed  residuals  of  the 
DBHO  equation  (Evans  et  al.  1967).  An  approx- 
imation to  the  normal  distribution  function 
may  also  be  used  (Burr  1967). 

Rethinnings  increase  d.b.h.  an  average  of 
0.4  inch  when  stand  densities  approximate 
growing  stock  levels  most  likely  to  be  goals 
of  timber  management. 


The  program  can  be  adapted  readily  for 
simulations  with  species  other  than  the  one 
used  for  the  test  problem.  Replacement  of 
tables  1  to  4,  several  statements  of  YIELD,  and 
species  designations  in  two  table  headings  are 
all  the  substitutions  necessary  if  the  same 
projection  procedure  is  used.  Statements  to 
be  replaced  are  the  equations  for  TOTO,  TOTT, 
DBHP,  PDBHE,  and  DBHO.  Projections  based 
on  other  procedures,  such  as  direct  estimation 
of  volume,  can  be  written  as  a  new  subroutine 
YIELD.  It  is  necessary  only  to  transfer  values 
of  CFMO(I),  BDFO(I),  CFMC(I),  and  BDFC(I) 
to  the  next  subroutine,  and  to  make  appropriate 
changes  in  the  READ  statements  in  INPUT1. 

Subroutine  ANVOL 

Subroutine  ANVOL  is  called  once  each  test 
to  compute  volume  per  acre  for  each  year  from 
initial  thinning  to  maximum  stand  age.  Volumes 
in  cubic  feet  and  board  feet  are  computed  by 
linear  interpolation  and  printed  in  a  composite 
table.  Stand  ages  cannot  exceed  179  years 
unless  dimensions  of  the  180-location  arrays  of 
acres  and  annual  volumes  are  increased. 

The  last  few  statements  of  the  routine 
expand  the  array  of  volumes  cut  to  assign  the 
volume  of  each  intermediate  cut  to  each  of  the 
years  before  the  next  cut  is  made.  These 
amounts  are  added  to  potential  postthinning 
volumes  in  subroutine  YEARS  to  compute 
volumes  per  acre  of  any  stands  not  given 
intermediate  cutting  because  they  are  older  than 
minimum  age  for  harvest. 

Subroutine  INPUT2 

INPUT2  is  called  once  each  game  to  enter 
numerical  values  of  variables  that  may  differ 
for  each  game  of  a  test.  Descriptive  data 
include  area  of  the  working  circle,  distribution 
of  area  by  age  classes,  nonstocked  area,  and 
number  of  acres  to  be  planted  annually. 
Various  costs  and  the  rate  at  which  they 
change  from  year  to  year  are  entered  into 
computer  memory.  One  to  ten  combinations 
of  price  limit,  allowable  cut,  and  minimum 
cutting  age  are  read  in  for  the  determination 
of  annual  cut  described  in  the  section  on  the 
data  deck. 

Values  read  by  INPUT2  that  do  not  appear 
in  other  tables  are  printed  as  a  record  of  game 
conditions  on  pages  headed  "alternatives  for 
this  game." 
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Subroutine  AREAS 

AREAS  is  called  once  each  game  to  compute 
volumes  and  area  distributions  at  the  end  of 
the  year  before  simulation  begins.  Acres  in 
each  1-year  age  class  are  expanded  to  obtain 
a  record  of  the  age  of  each  individual  acre. 
Total  area  (LAND)  cannot  exceed  1,000  acres 
unless  the  dimension  of  ACAGE(I)  is  increased. 
Age  of  the  oldest  acre  cannot  exceed  179  years 
unless  dimension  changes  listed  in  the  descrip- 
tion of  ANVOL  are  made.  Acres  are  totaled 
by  10-year  age  classes.  A  table  of  initial  dis- 
tribution of  acres  by  1-year  and  10-year  age 
classes  is  printed. 

Growing  stock  volume  is  totaled  in  board 
feet  and  in  cubic  feet.  Volume  of  an  acre  will 
be  added  to  the  total  of  only  one  of  the  two 
volume  units.  The  unit  will  be  board  feet  if 
the  board-foot  volume  equals  or  exceeds  the 
value  of  variable  BFMRCH  read  by  INPUT1. 
No  volume  will  be  computed  for  the  acre  if 
stand  age  is  less  than  the  specified  minimum 
(AGMRCH). 

Values  of  volume  and  money  variables 
that  are  not  zero  at  start  of  simulation  are 
computed.  These  values  are  then  assigned 
space  in  one  of  two  2-dimensional  arrays  to 
preserve  them  for  printing  by  OUTPUT2. 


Subroutine  YEARS 

YEARS  is  called  each  year  of  each  game  to 
simulate  the  changes  in  volume  and  value 
produced  by  tree  growth,  cutting,  and  other 
events.  Changes  and  their  order  of  occurrence 
are  indicated  by  the  comment  statements  of 
the  source  program.  Several  items  are  described 
in  more  detail  below. 

A  specified  area  (IPLNT)  is  planted  each 
year  if  nonstocked  acres  exist.  Nonstocked 
acres  are  those  deforested  by  fire  or  other 
catastrophe,  and  do  not  include  harvested  acres 
that  will  restock  in  the  allotted  time.  Some 
or  all  harvested  acres  could  be  added  to  non- 
stocked  area  to  simulate  delays  or  failures  in 
natural  regeneration. 

Age  of  each  acre  destroyed  and  added  to 
nonstocked  area  is  selected  at  random  with  a 
pseudorandom  number  generator  of  the  form: 


xi  ~  Ax-j_!  +  c  (modulo  M) 


(Greenberger  1961).  All  values  are  preset  except 
for  X  j  .  j  which  is  designated  as  variable 
ANUL.  The  generator  has  a  periodicity  of 
128.  Values  of  ANUL  from  0  to  127  may  be 
selected  at  random  to  vary  the  pattern  of  loss. 

Annual  harvest  in  acres  equals  the  constant 
or  variable  allowable  cut  less  any  losses  of 
whole  acres.  Volume  and  value  of  shelterwood 
or  seed  trees,  if  any,  are  credited  to  the  year  of 
final  cut  FINL  years  after  the  regeneration 
cycle  starts.  The  volume  may  increase  or 
decrease  during  the  regeneration  period,  and 
may  be  left  unharvested.  Desired  results  are 
obtained  by  entry  of  appropriate  values  for 
GROW  and  SHELT. 

Subroutine  YEARS  was  written  to  contain 
a  series  of  dynamic  events  useful  for  many 
species  and  forest  regions.  The  computations 
can  easily  be  added  to  or  modified  to  meet 
local  needs  or  to  test  special  alternatives.  Un- 
wanted alternatives  in  the  program  need  not  be 
removed.  They  can  be  bypassed  by  entry  of 
appropriate  values  for  variables  not  needed. 


Subroutine  OUTPUT1 

Subroutine  OUTPUT1  is  called  after  the 
first  year  of  each  game  and  at  the  end  of  each 
decade.  Numbers  of  acres  by  1-year  and  10- 
year  stand  age  classes  are  printed.  The  tables 
correspond  to  that  printed  by  AREAS  just 
prior  to  start  of  a  game. 


Subroutine  ANUAL 

ANUAL  is  called  every  year  of  each  game 
to  compute  40  volume,  area,  or  money  totals 
and  to  store  them  for  later  use.  Each  total 
is  stored  in  one  of  two  2-dimensional  arrays. 
The  first  dimension  identifies  the  variable,  the 
second  the  year  of  a  game  to  which  the  value 
applies.  Numerical  value  of  each  year  subscript 
is  year  plus  one  so  year  zero  of  a  game  can 
be  included  in  the  array.  Array  values  are 
used  in  all  subsequent  subroutines  of  MANAGD. 


Subroutine  OUTPUT2 

OUTPUT2  is  called  at  the  end  of  each  game 
to  print  the  results  of  each  year  of  the  game. 
Array  values  computed  and  stored  by  ANUAL 
are  printed  in  40  numbered  columns  that  extend 
across  four  pages  of  standard  Z-fold  paper. 
Entries  under  column  headings  are  printed  at 
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the  rate  of  40  lines,  or  years,  per  page.  Column 
headings  on  the  pages  produced  by  the  test 
problem  (appendix  2)  and  the  variable  lists  in 
the  source  program  of  ANUAL  (appendix  1) 
identify  the  variables  reported. 


Subroutine  WORTH 

Subroutine  WORTH  is  called  at  the  end  of 
each  game  to  discount  all  costs  incurred  and 
all  income  received.  Value  of  the  growing 
stock  at  the  end  of  the  simulation  period  is 
discounted  to  beginning  of  the  period.  The 
program  discounts  each  future  value  at  each 
of  20  compound  interest  rates.  Rates  range 
from  1.0  to  10.5  percent  at  intervals  of  0.5 
percent.  The  limits  and  interval  can  be  changed 
by  replacement  of  statements  for  CRATE(I) 
and  CRATE(K)  near  the  beginning  of  the 
subroutine. 

WORTH  prints  a  table  that  gives  the  present 
value  of  each  of  the  following  for  each  discount 
rate:  (1)  future  growing  stock,  (2)  all  incomes, 
(3)  sum  of  items  1  and  2,  (4)  all  costs,  and 
(5)  item  3  minus  the  sum  of  item  4  and  the 
value  of  the  growing  stock  at  beginning  of  the 
game.  Net  discounted  revenues  (present 
worths,  item  5)  may  be  plotted  over  discount 
rates  to  determine  the  internal  rate  of  return 
applicable  to  the  duration  and  conditions  of  the 
game. 


Subroutine  SUMRY 

Subroutine  SUMRY  may  be  called  at  the 
end  of  each  test  to  summarize  results  of  the 
games  of  the  test.  If  this  option  is  used, 
SUMRY  is  also  called  at  the  end  of  each  game 
to  store  specified  volume  and  money  values  in 
a  3-dimensional  array.  Values  stored  correspond 
to  the  columns  of  OUTPUT2  that  have  the 
column  numbers  entered  as  KOL(I)  by  INPUT1. 
Any  of  the  40  numbered  columns  of  OUTPUT2 
(appendix  2)  may  be  reproduced.  Not  more 
than  six  columns  may  be  summarized  during  one 
test  unless  the  dimensions  of  variables  KOL(I) 
and  SUMM(I,J,K)  are  increased.  Results  of  as 
many  as  10  games  may  be  summarized  at  one 
time. 

Results  of  the  games  of  a  test  are  printed 
together,  with  a  separate  page  for  each  variable 
selected  in  advance. 


Data  Deck 

Twenty-two  types  of  punch  cards,  listed 
below,  are  used  to  enter  initial  values  of  vari- 
ables into  computer  memory.  Most  cards  are 
not  optional  and  must  be  included  in  the  data 
deck  so  READ  statements  will  be  executed 
properly.  Four  types  are  optional  (6,  9, 10,  17) 
and  are  omitted  from  the  data  deck  if  the 
options  are  not  to  be  exercised. 

Data  cards  are  read  by  four  routines  in  the 
order  in  which  the  types  are  numbered.  Two 
types  are  read  once  by  the  main  program:  (1) 
card  type  1  or  BATCH  name,  and  (2)  card  type 
2  or  the  number  of  tests  to  be  performed  in  a 
batch.  These  identify  the  job  and  control  the 
number  of  times  the  rest  of  the  main  routine 
is  repeated. 

Card  types  3  to  10,  inclusive,  are  read  by 
INPUT1.  One  card  of  each  type  except  types 
6,  9,  and  10  must  be  read  once  each  test.  Card 
types  6  (1  card),  9  (15  cards),  and  10  (15  cards) 
are  omitted  from  the  data  deck  if  their  options 
are  not  to  be  used.  The  READ  statement  for 
card  type  6,  column  numbers  of  items  to  be 
summarized  by  SUMRY,  is  bypassed  when  zero 
is  punched  for  NKOLS  on  card  type  4.  Non- 
zero stumpage  prices  (BDPRI  and/ or  CFPRI) 
on  card  type  8  cause  the  corresponding  READ 
statements  for  variable  prices  of  card  types 
9  or  10  to  be  skipped. 

Subroutine  YIELD  reads  four  card  types 
once  each  test.  Types  11  to  14,  inclusive, 
contain  the  values  of  tables  1  to  4  of  this 
publication,  or  equivalent  information  applicable 
to  other  species  or  utilization  standards.  A 
type  consists  of  three  to  eight  punch  cards. 
Change  in  the  lowest  site  index  of  card  type 
12  (table  2)  or  in  the  minimum  diameter  of 
types  13  or  14  (tables  3  and  4)  will  require 
changes  in  the  statements  that  compute  array 
subscripts  and  probably  in  dimensions  of  the 
arrays. 

Card  types  15  to  22,  inclusive,  are  read  by 
INPUT2  once  each  game.  Each  type  consists 
of  one  card  except  for  optional  type  17,  which 
requires  10  punch  cards.  Statements  that  refer 
to  card  type  17,  variable  area  by  age  classes, 
are  bypassed  when  a  non-zero  value  is  punched 
for  KAREA  on  card  type  16. 

Card  types  18,  19,  and  20  contain  values  for 
the  price  control  procedure  of  Gould  and  O'Regan 
(1965).  The  number  of  acres  harvested  annually 
and  the  minimum  cutting  age  can  be  made  to 
vary  with  the  current  stumpage  price  of  saw 


-  8  - 


logs.  For  example,  in  the  second  game  of  the 
test  problem  (appendix  2),  5  acres  will  be  cut 
if  price  per  thousand  board  feet  does  not 
exceed  $12.00.  Seven  acres  will  be  cut  if  stump- 
age  price  is  $12.01  to  $15.00,  and  10  acres  will  be 
cut  if  price  exceeds  $15.00  but  is  less  than 
$99.00.  Sequence  of  harvest  is  from  oldest 
acre  to  youngest,  so  full  allowable  cut  will  be 


taken  only  if  sufficient  acres  above  minimum 
cutting  age  are  available.  If  price  control  is 
not  wanted,  entries  for  allowable  cut  in  columns 
1  to  4  of  card  type  19  and  for  cutting  age  in 
columns  1  to  8  of  card  type  20  are  the  desired 
constant  limits.  A  critical  price  greater  than 
the  largest  possible  price  (for  example,  $99.00) 
is  entered  in  columns  1  to  8  of  card  type  18. 


Order  and  Contents  of  the  Data  Deck 


Card  Read      Frequency    No.  of  Variable 

type    Optional        by  read  cards  name  Columns     Format  Description  of  variable 


1  NO  Main  Batch  1          BATCH (I)  1-24      3A8  Descriptive  name  to  identify  output  of 

one  pass  through  the  computer. 

2  NO  Main  Batch  1  NTSTS  1-4        14  Number  of  tests  in  the  batch,   each  with 

a  yield  table. 

3  NO  INPUTl        Test  1  DESCR(I)  1-40      5A8  Phrase  to  describe  conditions  of  one 

test;   to  identify  output. 

4  NO  INPUTl        Test  1          NGAME  1-4        14  Number  of  trials  (games)  to  be  operated 

in  one  test. 

NKOLS  5-8        14  Number  of  columns  of  OUTPUT2  to  be 

printed  by  SUMRY. 

NOYRS  9-12      14  Number  of  years  simulated  in  each  game. 

Can  be  up  to  150,  but  will  usually  be 
less . 


INPUTl        Test  1         AGEO  1-8        F8.3        Stand  age  at  time  of  initial  thinning. 

First  age  given  in  the  yield  table. 

SITE  9-16      F8o3        Site  index.     Base  age  and  crown  class 

same  as  used  to  derive  growth  equations. 

DENO  17-24      F8.3        Number  of  trees  per  acre  at  age  AGEO. 

DBHO  25-32      F8.3        Average  diameter  breast  high  of  the 

stand  at  age  AGEO. 

ROTA  33-40      F8.3        Maximum  age  in  the  yield  table;   1  year 

more  than  the  maximum  age  expected  dur- 
ing simulations.     Cannot  exceed  180 
years . 

PRET  41-48      F8.3        Estimated  percentage  of  the  number  of 

trees  to  be  retained  in  initial  thinning 
at  age  AGEO.     Enter  as  a  percent,  e.g., 
35.0. 

DLEV  49-56      F8.3        Density  level  for  intermediate  cuts 

after  initial  thinning.     Based  on  table 
1  of  this  publication  and  procedure 
described  in  YIELD. 

CYCL  57-64      F8.3        Interval  between  intermediate  cuts. 

Equal  to  or  a  multiple  of  RINT. 

RINT  65-72      F8.3        Number  of  years  for  which  a  growth  pro- 

jection is  made  by  the  equation  in  YIELD. 

THIN  73-80      F8.3        Density  level  after  initial  thinning  at 

age  AGEO.     Based  on  table  1  and  proce- 
dure described  in  YIELD.     May  equal  DLEV. 
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Card  Read      Frequency    No.  of  Variable 

type    Optional        by  read  cards  name  Columns  Format 


Description  of  variable 


6 

YES 

INPUT 1 

Test 

1 

KOL(l) 

1-4 

14 

K0L(2) 

5-8 

14 

KOL(3) 

9-12 

14 

KOL(4) 

13-16 

14 

K0L(5) 

17-20 

14 

KOL(6) 

21-24 

14 

7 

NO 

INPUT 1 

Test 

1 

AGMRCH 

1-8 

F8.3 

BFMRCH 

9-16 

F8.3 

17-24 

F8 . 3 

to.  J 

COMBF 

33-40 

F8 . 3 

BFPCT 

41-48 

F8.3 

CFPCT 

49-56 

F8.3 

GNTR 

57-64 

F8.3 

8 

NO 

INPUT 1 

Test 

1 

BDPRI 

1-8 

F8.3 

CFPRI 

Q    1  A 

F8 . 3 

9 

YES 

INPUT  1 

Test 

15 

PRICF(I) 

1-80 

10F8, 

10 

YES 

INPUT  1 

Test 

15 

PRIBD(I) 

1-80 

10F8. 

11 

NO 

YIELD 

Test 

4 

TABLl(K) 

1-63 

21F3. 

12 

NO 

YIELD 

Test 

3 

TABL2(K,L) 

1-75 

25F3, 

Numbers  of  the  columns  of  OUTPUT 2  to  be 
printed  by  SUMRY.     Maximum  number  of 
columns  is  6  and  must  agree  with  NKOLS 
of  card  type  4.     Column  numbers  are  from 
1  to  40,  as  given  in  the  column  headings 
of  the  printout  of  0UTPUT2  of  the  test 
problem. 

Minimum  stand  age  for  an  acre  to  be  in- 
cluded in  growing  stock  volume. 

Minimum  volume  in  M  bd.  ft.  for  an  acre 
to  be  included  in  board-foot  growing 
stock  volume. 

Minimum  volume  per  acre  in  M  bd.  ft.  for 
commercial  salvage  after  fire,  wind,  or 
other  loss. 

Minimum  cut  per  acre  in  merchantable 
cubic  feet  (table  3)  for  a  cut  to  be  of 
positive  commercial  value. 

Minimum  cut  per  acre  in  M  bd.  ft. 
(table  4)  for  a  cut  to  be  of  positive 
commercial  value. 

Ratio,  as  a  decimal,  of  board-foot 
stumpage  values  of  thinnings  to  board- 
foot  stumpage  values  of  harvests. 

Ratio,  as  a  decimal,  of  cubic-foot 
stumpage  values  of  thinnings  to  cubic- 
foot  stumpage  values  of  harvests. 

Any  number  between  0  and  1023  used  to 
generate  random  element  of  the  increase 
from  DBHT  to  DBH0.     Enter  number  larger 
than  1024  to  bypass  this  step. 

Stumpage  price  per  M  bd.  ft.  of  final 
harvest  if  price  is  constant  for  all 
years  of  a  game.     Enter  zero  if  variable 
prices  will  be  entered  with  card  type  10. 

Stumpage  price  per  100  cubic  feet  of 
final  harvest  if  price  is  constant  for 
all  years  of  a  game.     Enter  zero  if  vari- 
able prices  will  be  entered  with  card 
type  9. 


Stumpage  price  per  100  cubic  feet  of 
harvest  for  each  of  150  years.  Used 
when  CFPRI  equals  zero. 


Stumpage  price  per  M  bd, 
for  each  of  150  years, 
equals  zero. 


ft.  of  harvest 
Used  when  BDPRI 


1  Basal  area  after  thinning  in  relation  to 
stand  diameter.  Values  of  table  1  copied 
on  punch  cards.     Used  with  DLEV  and  THIN. 


Tree  heights  by  age  and  site  index  class. 
Values  of  table  2  copied  on  punch  cards. 
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Card  Read    Frequency    No.   of  Variable 

type    Optional  by  read  cards  name  Columns  Format 


Description  of  variable 


13  NO  YIELD          Test  8  TABL3 (K)  1-72      24F3.3    Factors  for  conversion  of  total  cubic 

feet  to  merchantable  cubic  feet.  Values 
of  table  3  copied  on  punch  cards. 


14          NO  YIELD         Test  6  TABL4 (K)  1-78      26F3.2    Factors  for  conversion  of  total  cubic 

feet  to  board  feet.  Values  of  table  4 
copied  on  punch  cards. 


15  NO  INPUT2        Game  1          GMNAM(I)  1-24      3A8  Descriptive  name  to  identify  each  game 

of  a  test. 

16  NO  INPUT 2        Game  1          LAND  1-4        14  Total  acres  in  simulated  working  circle. 

Maximum  is  1,000  acres. 


MOLD  5-8        14  Age  of  oldest  stand  in  the  working  cir- 

cle at  start  of  a  game.  Maximum  is  179 
years . 

NONSTK  9-12      14  Number  of  acres  non-stocked  at  start  of 

a  game.     Does  not  include  acres  harvest- 
ed the  year  before  simulation  begins  if 
regeneration  will  take  place  in  the 
allotted  time. 


KAREA  13-16      14  Number  of  acres  in  each  1-year  age  class 

when  there  is  equal  area  in  each  class 
except  for  NONSTK. 


17-20      14  Number  of  acres  of  NONSTK  regenerated 

annually  by  direct  seeding  or  planting 
at  a  cost  of  CPLT  per  acre. 


17          YES  INPUT2        Game  10  IACRE(I)  1-72       1814        Acres  in  each  1-year  age  class  from  0  to 

not  more  than  179.    Use  if  constant  area 
KAREA  is  not  wanted.     Include  NONSTK  in 
IACRE(l)  as  well  as  on  card  type  16. 


18  NO  INPUT2        Game  1  PRIDIV(I) 

19  NO  INPUT2        Game  1  MALCUT(I) 

20  NO  INPUT2        Game  1  FMRCHD(I) 


1-80      10F8.3    Limiting  prices  used  to  determine  annual 
cut  in  acres  and  minimum  cutting  age. 

1-40      1014        Allowable  annual  cut  in  acres.  May 
vary  with  PRIDIV. 

1-80      10F8.3    Minimum  cutting  age.     May  vary  with 
PRIDIV. 


21         NO  INPUT2        Game  1  SHELT  1-8        F8.3        Volume  of  shelterwood  in  M  bd.  ft.  Enter 

zero  if  shelterwood  or  seed  trees  are 
not  retained. 

RATE  9-16      F8.3        Rate  of  annual  increase  in  costs.  Enter 

zero  if  constant  costs  are  desired. 
Otherwise,  enter  percentage  as  a  decimal. 

CPLT  17-24      F8.3        Cost  of  regenerating  1  acre  by  seeding 

or  planting. 

CTHN  25-32      F8.3        Cost  per  acre  of  precommercial  thinning 

with  stand  conditions  as  specified  for 
the  simulation. 

CLOSS  33-40      F8.3        Cost  per  acre  of  cleanup  after  loss  due 

to  fire,  wind,  etc.,  when  volume  that 
can  be  salvaged  is  less  than  BFSALV. 
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Card  Read      Frequency    No.  of  Variable 

type    Optional        by  read  cards  name  Columns  Format 


Description  of  variable 


ACCST 


41-48  F8.3 


Total  per  acre  for  1  year  of  the  annual 
costs  that  can  be  assessed  by  area. 


CUCST 


49-56  F8.3 


Total  of  the  costs  that  can  be  assessed 
against  each  100  cubic  feet  harvested. 


57-64  F8.3 


Total  of  the  costs  that  can  be  assessed 
against  each  M  bd.  ft.  harvested. 


65-72      F8.3        Percentage  increase  or  decrease  in 

shelterwood  volume  during  the  regenera- 
tion period.     Enter  as  a  decimal. 
Enter  -1.0  if  shelterwood  or  seed  trees 
will  not  be  harvested. 


73-80      F8.3        Number  of  years  between  harvest  cut  and 
removal  of  shelterwood  or  seed  trees. 
Enter  zero  if  not  to  be  removed. 


22 


INPUT2  Game 


F8.5 


Percentage,  as  a  decimal,  of  the  area  of 
forest  lost  annually  to  fire,  wind,  etc. 


9-16      F8.5        Any  number  between  0  and  127  used  to 

begin  generation  of  pseudorandom  numbers 
that  represent  ages  of  stands  lost  to 
fire  or  other  agency. 


Description  of  Test  Problem 

The  test  problem  that  follows,  (detailed  in 
appendix  2)  demonstrates  most  computations 
possible  and  the  printed  results  obtained.  It 
may  be  used  to  verify  accuracy  of  punching  of 
source  decks  and  compatibility  of  the  program 
with  locally  available  compilers.  Growth  projec- 
tions and  volume  conversions  are  based  on 
relationships  applicable  to  Black  Hills  ponderosa 
pine  (Myers  1966).  Data  for  costs,  stumpage 
prices,  and  other  items  are  hypothetical.  Results 
of  the  simulations  are  therefore  examples  only, 
and  do  not  apply  to  any  real  forest  area. 

Assume  an  area  of  915  acres  of  managed 
stands  that  range  from  0  (just  harvested)  to 
129  years  old.  There  are  7  acres  in  each  1-year 
age  class,  plus  5  acres  of  old  burn  and  wind- 
throw  that  must  be  seeded  or  planted.  Annual 
losses  to  fire,  wind,  and  other  agencies  average 
0.04  percent  of  the  forested  area.  Site  index 
of  all  acres  is  60  feet. 

Stands  will  be  regenerated  by  two-cut 
shelterwood,  and  will  be  thinned  several  times 
during  a  rotation.  Shelterwood  volume  will 
average  4,000  board  feet  per  acre,  and  will 


increase  an  average  of  3  percent  of  initial 
volume  each  of  the  10  years  before  the  final 
cut.  Stands  will  be  thinned  for  the  first  time 
when  they  are  30  years  old.  At  this  age,  there 
will  be  1,000  trees  per  acre  that  average  4.5 
inches  in  diameter.  Initial  thinning  will  be  to 
level  120,  or  120/80  times  the  basal  areas  in 
table  1.  Rethinnings  at  20-year  intervals  will 
be  to  level  100,  or  100/80  times  tabulated  basal 
areas.  Stands  are  not  expected  to  ever  become 
150  years  old. 

Potential  prices  of  two  products  have  been 
estimated  for  each  of  the  next  30  years.  The 
stumpage  price  of  100  cubic  feet  of  roundwood 
from  mature  trees  or  from  thinnings  is  expected 
to  be  $2.50  throughout  the  period.  Price  of  a 
thousand  board  feet  of  saw  logs  is  expected  to 
vary  annually,  as  shown  in  column  28  of  the 
printout  of  annual  results  (appendix  2).  Saw 
logs  from  thinnings  will  sell  for  85  percent  of 
the  price  of  logs  from  harvest  cuts.  A 
minimum  commercial  cut  of  saw  logs  will  be 
3,000  board  feet  per  acre,  except  that  1,500  board 
feet  may  be  salvaged  from  an  acre  after  fire 
or  other  catastrophe.  Minimum  commercial 
cut  of  roundwood  will  be  400  cubic  feet. 
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Current  value  of  the  growing  stock  will  be 
computed  only  for  stands  at  least  40  years  old. 
Value  will  be  computed  for  cubic  volume  for 
acres  with  less  than  1,500  board  feet.  Other- 
wise, board-foot  volumes  will  be  used. 

Present  costs  of  various  operations  are  as 
follows: 

Costs  per  acre- 
Seeding— $30 

Precommercial  thinning— $25 
Cleanup  where  salvage  is  not  pos- 
sible—$25 
Annual  costs— $0.20 
Costs  assessed  against  volume  sold— 
Per  100  cubic  feet— $0.05 
Per  thousand  board  feet— $1.56 

These  costs  are  expected  to  increase  at  a  rate 
of  1  percent  annually.  Resources  are  available 
to  seed  1  acre  each  year. 

Two  possible  means  of  setting  the  allowable 
annual  cut  are  to  be  tested.  One  alternative 
is  to  harvest  7  acres  annually,  regardless  of 
price  fluctuations.  Stands  less  than  130  years 
old  will  not  be  cut.  A  second  possibility  is  to 
harvest:  (1)  5  acres  if  stumpage  price  per 
thousand  board  feet  is  $12  or  less,  (2)  7  acres 
if  the  price  is  $12.01  to  $15,  and  (3)  10  acres 
if  price  exceeds  $15  per  thousand.  Stands  less 
than  130  years  old  will  not  be  cut  except  that 
the  minimum  age  will  be  120  years  when 
stumpage  price  exceeds  $15.  Periodic  production 
in  board  feet  and  net  income  will  be  compared. 
Values  needed  to  obtain  present  worths  will 
be  computed. 

Data  cards  to  enter  these  values  into  compu- 
ter memory  must  contain  the  alphameric 
characters  given  in  the  following  list.  Decimal 
points  are  shown  for  numbers  in  F-format  to 
indicate  the  way  in  which  percentages  and 
money  values  are  entered.  Spaces  between 
numbers  do  not  correspond  to  the  blank  columns 
of  the  data  cards.  Card  types  9  and  17  are  not 
included  in  the  data  deck  because  the  options 
that  require  them  will  be  bypassed. 

Test  conditions  and  results  of  the  simula- 
tions are  printed  on  seven  types  of  pages 
(appendix  2).  The  first  two  types,  (1)  a  yield 
table  and  (2)  tables  of  volumes  per  acre  for 
each  year  of  stand  age,  appear  once  because 
one  test  was  run.  Four  types  of  pages  are 
printed  for  each  of  the  two  games.  The 
seventh  type  of  page  appears  once  at  the  end 
of  the  printout  to  summarize  specified  results 
of  the  two  games. 


The  two  sheets  of  "alternatives  for  this 
game"  show  the  values  used  in  the  simulations, 
including  the  different  allowable  cuts  and  cutting 
ages  tested. 

Distributions  of  acres  by  age  classes  (page 
type  four)  appear  on  two  sets  of  pages,  one  set 
for  each  game.  Pages  for  year  zero  show  7 
acres  in  each  of  129  1-year  age  classes.  Age 
class  zero  has  an  additional  5  acres  of  non- 
stocked  area.  Acreages  are  the  same  for  both 
games,  because  initial  distributions  were  the 
same.  Similar  pages  are  printed  at  the  end  of 
the  first  year  of  each  game  and  at  the  end  of 
each  decade.  For  brevity,  only  the  page  printed 
after  the  thirtieth  year  of  each  game  is  repro- 
duced in  appendix  2.  After  30  years  of  simula- 
tion, losses  and  direct  seeding  have  modified 
the  pattern  of  7-acre  units.  In  addition,  area 
distribution  of  the  second  game  has  been 
changed  by  the  variable  annual  cuts. 

The  fifth  type  of  page  is  a  set  of  four  pages 
for  each  game.  Values  in  many  of  the  40 
numbered  columns  differ  between  games. 
Board-foot  volumes  are  unequal  because  of 
variations  in  annual  cuts  of  mature  timber 
during  the  second  game.  This  caused  money 
values  to  differ  from  those  reported  for  the 
first  game. 

A  page  of  discounted  money  values,  the 
sixth  type  of  page,  is  printed  for  each  game. 
Rate  of  return  was  about  the  same  for  both 
games.  Both  operations  were  profitable.  In 
addition,  the  forest  would  probably  be  in  good 
condition  to  produce  other  products,  especially 
recreation. 

Last,  specified  values  from  each  game  were 
printed  together  for  convenience  in  interpreta- 
tion of  results.  Total  volume  in  board  feet  of 
all  cuts  plus  growing  stock  (column  10)  was 
higher  after  30  years  where  equal  areas  were 
cut  each  year.  Differences  were  never  great; 
variable  annual  cuts  produced  the  larger 
volume  after  20  years.  Total  net  worth  (column 
40)  was  greater  where  annual  cuts  varied  with 
price,  except  for  several  of  the  earlier  years. 

It  must  be  emphasized  that  results  of  these 
or  other  simulations  depend  on:  (1)  duration  of 
the  games,  (2)  values  entered  for  the  various 
variables,  (3)  assumptions  made,  and  (4)  degree 
to  which  the  model  represents  reality. 

The  above  information,  additional  data,  and 
knowledge  of  local  conditions  would  help  the 
forest  manager  decide  how  he  might  best  con- 
duct his  business.  Money  yields  would  en- 
courage the  manager  to  vary  annual  cuts  in 
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Data  Deck  for  Test  Problem 


Card 
type 

Contents 

of  Cards 

1 

TEST  PROBLEM 

2 

1 

3 

MANAGED , 

T  U  T  VTM  TP  T\ 

IhlNNEU 

AGE  30 

4 

2 

2 

30 

5 

30 . 

60. 

1000.  4.5 

150. 

40. 

100. 

20. 

10. 

120. 

6 

10 

40 

7 
8 

40 . 
U .  U 

1.5 
2.50 

1.5  400. 

3.0 

.85 

1. 

0 

2222. 

10 

(Column 

28  of  printout  by  0UTPUT2) 

11 

(Table  1 

of  this 

publication) 

12 

(Table  2 

of  this  publication) 

13 

(Table  3 

of  this 

publication) 

14 

(Table  4 

of  this  publication) 

15 

EQUAL  AREAS  CUT  ANNUALLY 

16 

915 

129 

5  7 

1 

18 

99. 

0 

0  0 

0 

0 

0 

0 

0 

0 

19 

7 

0 

0  0 

0 

0 

0 

0 

0 

0 

20 

130. 

0 

0  0 

0 

0 

0 

0 

0 

0 

21 

4.0 

.01 

30.  25. 

25. 

.20 

05 

1. 

56 

30 

10.0 

22 

.0004  21. 

15 

VARY  CUT 

WITH  PRICE 

16 

915 

129 

5  7 

1 

18 

12. 

15. 

99.  0 

0 

0 

0 

0 

0 

0 

1  Q 

5 

7 

10  0 

o 

o 

o 

o 

o 

o 

20 

130. 

130. 

120.  0 

0 

0 

0 

0 

0 

0 

21 

4.0 

.01 

30.  25. 

25. 

.20 

05 

1. 

56 

.30 

10.0 

22 

.0004  21. 
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response  to  changes  in  stump  age  price.  Highly 
variable  annual  cuts  and  equally  variable  net 
incomes  would  suggest  that  additional  simula- 
tions be  run  to  test  other  alternatives.  Cost 
of  computer  time  need  not  restrict  the  manager 
in  his  search  for  information.  The  test  problem 
was  compiled  and  run  in  26  seconds  of  central 
processor  time  and  9  seconds  of  input-output 
time. 
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APPENDIX  I 
Program  MANAGD 


DEFINITIONS   OF  VARIABLES 

AC  AG  E ( I  )    =   ONE   ACRE,    WITH   AGE   OF   NUMERICAL   VALUE    AND    IDENTIFIED  BY 

SUBSCRIPT 
ACCST  =   ANNUAL   COST   PER  ACRE 
AG'EO   =    INITIAL   AGE    IN   YIELD  TABLE 

AGMRCH   =   MINIMUM  AGE   FOR   STAND   TO   BE    INCLUDED    IN   GROWING  STOCK 
ANBDFI I)    =   M   BD.    FT.    PER   ACRE   AT   ENO   OF   EACH  YEAR 
ANCUV(I)    =   CU.FT.    STANDING   PER   ACRE   AT   END   OF   EACH  YEAR 
ANNET   =   ANNUAL   NET  INCOME 

ANUL   =   NUMBER   BETWEEN   0  AND   127   USED   TO   START  GENERATION  OF 

PSEUDORANDOM  NUMBERS 
BASC   =    BASAL   AREA   REMOVED   PER  ACRE 
BASO   =   BASAL    AREA   PER   ACRE   BEFORE  THINNING 
HAST   =   BASAL   AREA   PER   ACRE    AFTER  THINNING 
BATCH ( I )    =   JOB  NAME 

BDFC(I)    =   M   BD.    FT.    REMOVED   PER  ACRF 
BDFn ( I )    =   M   BD.    FT.    PER   ACRE    BEFORE  THINNING 
BDFT(I)    =   M   BD.    FT.    PER   ACRE   AFTER  THINNING 
BDPRI    =   CONSTANT   STUMPAGE    PRICE   PER   M   BD.  FT. 
BFCST   =   COSTS   PER   M   BD.    FT.  HARVESTED 
BFMRCH   =   MINIMUM   VOLUME   TO   BE    INCLUDED    IN   BD. FT , 


GROWING  STOCK 


BFPCT   =    PCT.    TO   CONVERT    BD.FT.    PRICE   FOR  THINNINGS 

BFSALV  =   MINIMUM   BD.FT.    FOR   COMMERCIAL  SALVAGE 

CFMC(I)    =   MERCHANTABLE   CU.FT.    REMOVED   PER  ACRE 

CFMO(I)    =   MERCHANTABLF   CU.FT.    PER   ACRE    BEFORE  THINNING 

CFMT(I)    =   MERCHANTABLE   CU.FT.    PER   ACRE    AFTER  THINNING 

CFPCT   =   PCT.    TO  CONVERT   CU.FT.    PRICE   FOR  THINNINGS 

CFPRI    =   CONSTANT   STUMPAGE    PRICE    PER   100  CU.FT. 

CLOSS   =   COST   OF   CLEANUP   OF   VOLUME   NOT  SALVAGED 

COMBF  =   MINIMUM   COMMERCIAL   CUT    IN   BD.  FT. 

COMCU  =    MINIMUM   COMMERCIAL   CUT    IN   CU.  FT. 

CPLT   =    PLANTING     COST   PER  ACRE 

CRATE (  I  )   =   INTEREST  RATES   FOR  DISCOUNTING 

CSTAC   =   ANNUAL   COSTS   BASED  ON  AREA 

CSTVL   =    ANNUAL   COSTS   FOR   VOLUME  HARVESTED 

CTHN   =   COST    PER   ACRE   OF    PR ECOMM ERC I AL  THINNING 

CUCST  =  COSTS  PER   100  CUBIC   FEET  HARVESTED 

CUTAGE   =   MINIMUM  CUTTING  AGE 

CYCL   =    INTERVAL   BETWEEN    INTERMEDIATE  CUTS 

QRHO   =    AVERAGE    STAND   D.B.H.    BEFORE  THINNING 

DBHT   =   AVERAGE    STAND   D.B.H.    AFTER  THINNING 

DEFOR   =    PERCENTAGE ( AS   DECIMAL)    OF   NUMBER   OF   ACRES   LOST  ANNUALLY 

DENC   =   TREES   REMOVED   PER  ACRE 

DENG  =  TREES   PER  ACRE   BEFORE  THINNING 

DENT   =   TREES   PER   ACRE   AFTER  THINNING 

DESCR(I)    =   DESCRIPTION  OF  TEST  CONDITIONS 

DISC!  I)    =  DISCOUNTED  VALUE   OF   FUTURE  COSTS 

DISG(I)    =   DISCOUNTED   VALUE    OF   GROWING  STOCK 

DISI(I)    =   DISCOUNTED  VALUE   OF   FUTURE  INCOMES 

DLEV   =   GROWING   STOCK    LEVEL    FOR   SECOND   AND   SUBSEQUENT  THINNINGS 
FINL   =   YEARS   BETWEEN   HARVEST   AND  REMOVAL   OF  SHELTERWOOD 
FMRCHD ( I )   =  MINIMUM  CUTTING  AGE   BASED  ON  PRICE 
GMNAM(I)    =   NAME   OF   THE  GAME 

GNTR  =    PSEUDORANDOM   NUMBER   GENERATOR.    VALUE    0   TO  1023. 

GROW   =   GROWTH   RATE   OF    SHELTER  WOOD 

GSVALB  =  DOLLAR  VALUE   OF   BD.FT.   GROWING  STOCK 

GSVALC  =  DOLLAR  VALUE   OF   CU.FT.   GROWING  STOCK 

GVLBF  =  GROWING  STOCK  VOLUME,    M  BD.  FT. 

GVLCU   =    GROWING    STOCK    VOLUME,  CU.FT. 

HTSO   =   TREE   HEIGHT   BEFORE  THINNING 

HTST   =   TREE   HEIGHT    AFTER  THINNING 

lACRE(I)    =   ACRES   OF   WORKING   CIRCLE   IN   EACH   1-YEAR    AGE   CLASS,  AT 

START   OF  GAME 
IALCUT  =   NUMBER   OF   ACRES   ALLOWABLE   ANNUAL  CUT 
IGAME   =   NUMBER   OF  GAME 

IPLNT   =   NUMBER   OF   NON-STOCKED   ACRES   REGENERATED  ANNUALLY 
ISUM(I)    =   TOTAL   ACRES   EACH   10-YR   AGE  CLASS 
ITEST   =   NUMBER   OF  TEST 

IVARU.J)    =   VARIABLES    TO   BE   PRINTED   BY  0UTPUT2 

IYEAR   =   YEAR   WITHIN   RUN   OF    A  GAME 

KAREA  =    EQUAL   AREA   OF    EACH    1-YEAR   AGE  CLASS 

KOL(I)    =   COLUMN   NUMBER ( FROM   0UTPUT2)    PRINTED   BY  SUMRY 

KOUNT   =   COUNT   OF   ACRES   HARVESTED,    PLUS  ONE 

LAND  =   TOTAL   ACRES    IN   SIMULATED   WORKING  CIRCLE 

LAST   =   NUMBER   OF   LAST   ACRE  HARVESTED 

MALCUT(I)    =   ANNUAL   ALLOWABLE   CUT   BASED   ON  PRICE 

MOLD  =  AGE  OF  OLDEST  ACRE  IN  WORKING  CIRCLE  AT  START  OF  A  GAME 
UACREII)    =   ACRES   OF   WORKING   CIRCLE    IN   EACH    1-YEAR   AGE  CLASS, 

DURING  A  GAME 
NGAMF   =   NUMBER   OF   GAMES   PER  TEST 

NKOLS   =    NUMBER   OF   COLUMNS   OF   0UTPUT2    TO   BE    PRINTED   BY  SUMRY 

NONSTK  =  NONSTOCKED  AREA 

NOYRS  =   NUMBER   OF  YEARS    IN   A  GAME 

NTSTS   =   NUMBER  OF   TESTS    IN  BATCH 

PRET  =   PERCENTAGE   OF   TREES   RETAINED  AFTER   INITIAL  THINNING 
PREVII)   =  PRESENT  VALUE   OF  GROWING  STOCK  AND  INCOMES 
PRIBD(I)    =   VARIABLE    STUMPAGE   PRICE   PER   M   BD.  FT. 
PRICF(I)    =   VARIABLE   STUMPAGE   PRICE   PER    100  CU.FT. 
PRIDIVII)   =   PRICES  USED  TO   SET  POLICY 
PWTH( I )    =   PRESENT  WORTH 

RATE   =   RATE   OF   ANNUAL    INCREASE    IN  COSTS 
RETHV   =   ANNUAL   RETURN   FROM   FINAL  HARVEST 
RETRN  =   ANNUAL    INCOME   FROM  STUMPAGE 
RETTH   =   ANNUAL   RETURN   FROM  THINNINGS 

RINT   =   NUMBER   OF   YEARS   FOR   WHICH   GROWTH   PROJECTION    IS  MADE 

ROTA  =  LONGEST  POSSIBLE   ROTATION    IN  YIELD  TABLE 

SCPLT   =   TOTAL   ANNUAL   PLANTING  COST 

SCTHN  =   SUM  OF   PRECOMMERC I AL   THINNING  COSTS 

SHELT   =   M   BD.    FT.    RETAINED   AS   SEED   TREES   OR  SHFLTERWOOD 

SHWD   =   VOLUME   OF   SHELTERWOOD   AT    FINAL  CUT 

SITE   =   SITE  INDEX 

SUMM(I,J,K)    =   ARRAY   FOR   PRINTING   BY  SUMRY 


TCOST   =   TOTAL   ANNUAL  COSTS 

THIN   =   GROWING   STOCK   LEVEL    FOR    INITIAL  THINNING 

TOTC    =   TOTAL   CUBIC    FEET   REMOVED   PER  ACRE 

TOTO  =  TOTAL   CUBIC   FEET   PER   ACRE    BEFORE  THINNING 

TOT T   =   TOTAL   CUBIC   FEET   PER   ACRE   AFTER  THINNING 

VARII.J)    =  VARIABLES   TO   BE   PRINTED   BY  0UTPUT2 

VBHV   =   BOARD-FOOT   VOLUME   FROM  HARVESTS 

VBTH   =   BOARD-FOOT   VOLUME   FROM  THINNING 

VCHV   =   CUBIC-FOOT   VOLUME   FROM  HARVESTS 

VCTH   =   CUBIC-FOOT   VOLUME   FROM  THINNING 

VLBF   =   VOLUME   HARVESTED,    M   BD.  FT. 

VLCU   =   VOLUME   HARVESTED,  CU.FT. 

YRLOS   =   NUMBER   OF   ACRES   LOST  ANNUALLY 

COMMON   ACCST, AGEO, AGMRCH, ANBDFI 181) , ANCUV! 181) , ANNET , BATCH ( 3 ), BDFC 
1(180) , BDFOl 180) , BFCST, CFMC (180) ,CFMO( 1 80 ) , C STAC , C ST VL , C UC ST , CUT AGE 
2,DBH0,DENO,DESCR( 5) , DL EV , FMRCHD ( 1 0 ) , GMN AM ( 3 ) , GSVALB , GSVALC , GVL BF , 
3G VLCU, I  ACRE! 180) , IALCUT, IGAME, I  SUM ( 18),ITEST,IVAR(26,150),IYEAR, 
4K0LI6) , LAND, LAST, MALCUTI 10) , MOLD , NACRE ( 1 80 ) , NG AME , NKOL S , Nl , NONS T K , 
5N0YRS,PRET,PRIBD( 150) ,PRICF(150),PRIDIV< 10  ), RETRN , ROTA , RATE , GROW, 
6SITE,SUMM(6, 25, 10 ), TCOST , SHELT , YRLOS , F I NL , CLOS S , CTHN , CPLT , C F PCT , 
7VAR(14,150)  ,VLBF, VLCU, AC AGE (  1000)  , C YCL , KOUNT , DEFOR , ANUL , M I  X , I  PL NT , 
8RINT , THIN, BFMRCH, BFSALV, COMCU, COMBF, BFPCT, GNTR 

READ    (5,1)    (BATCH!  I ), 1  =  1 ,3  I 

1  FORMAT  (3A8) 
READ    (5,2)  NTSTS 

2  FORMAT  (14) 

DO  8    ITEST=1, NTSTS 
CALL  INPUT1 

PRINT   YIELD   TABLE    AND   COMPUTE   VOLUME    FOR   EACH   YEAR   OF   STAND  AGE 

CALL  YIELD 
CALL  ANVOL 


OPERATE   SYSTEM    FOR   DESIRED  NUMBER   OF  GAMES 


CREATE    ACRES    IN   EACH   AGE  CLASS 


CALL  AREAS 
UK   =  0 


OPERATE  SYSTEM  FOR  DESIRED  NUMBER  OF  YEARS 
DO   7    I  YE  AR= 1 , NOYRS 


IF    I  I  YEAR    .LE.    1)    GO   TO  3 
IF    (UK   .EQ.    10)   GO  TO  4 

GO  TO  6 

3  UK  =  1 
GO  TO  5 

4  UK   =  0 

5  CALL  0UTPUT1 

6  UK   =    UK   +  1 
CALL  ANUAL 

7  CONTINUE 

PRINT   VOLUMES   AND  VALUES   FOR   EACH  YEAR 


SUMMARIZE   DESIRED   NUMBER   OF   COLUMNS   OF  0UTPUT2 

IF    (NKOLS    .LE.    0)    GO   TO  8 
CALL  SUMRY 
8  CUNTINUE 
CALL  EXIT 

END 

SUBROUTINE  INPUT1 

COMMON   ACCST, AGEO, AGMRCH, ANBDFI 181) , ANCUVI  181)  , ANNE T , BATCH ( 3 ) , BDFC 
1(180), BDF0I180) , BFCST, CFMC I  180) , C FMO I  1 80 ) , C ST AC , CST VL , CUC ST , CUT AGE 
2, DBH0,DEN0,DESCR(5), DLEV, FMRCHD! 10), GMNAM(3),GSVALB, GSVALC, GVLBF , 
3GVLCU, I  ACRE! 180) , IALCUT, I G AME , I  SUM ( 1 8 ) , I  TEST, IV AR( 26, 150), IYEAR, 
4K0L(6) , LAND, LAST, MALCUTI 10), MOLD, N ACRE!  1 80 ) , NGAME , NKOL S , N 1 , NDNSTK  , 
5 NOYRS, PRET, PRIBDI150) , PR ICF (150), PR  I DIV( 10), RETRN, ROTA, RATE, GROW, 
6SI TE, SUMM ( 6, 25, 10) , TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT, CFPCT, 
7VARI  14, 150), VLBF, VLCU, ACAGEI  1000)  , C YC L , KOUNT , DE F OR , ANUL , M I  X , I  PL NT  , 
8R I  NT, THIN, BFMRCH, BFSALV, COMCU, COMBF, BFPCT, GNTR 

SET    INITIAL   VALUES   OF  ZERO 

UO    1  1=1,6 

1  KOL(I)    =  0 
DO  2  1=1,150 
PR  I BD (  I  )    =  0.0 

2  PRICF ( I )    =  0.0 
DO   3  1=1,181 
ANBDFI I )    =  0.0 

3  ANCUVI I )    =  0.0 


DO  4  ! 
DO  4  . 
DO  4  t 
SUMM I i 


>25 


READ   VALUES   THAT   DO   NOT   CHANGE  WITHIN 


-  16  - 


READ   (5,5)  (DESCR(I),I=1,5) 

5  FORMAT  (5A8) 

READ    (5,6)    NGAME, NKOLS, NOYRS 

6  FORMAT  (2014) 

READ    (5,7)    AGEO, SITE, DENO.DBHO, ROTA, PRE T , DL E V, C YCL , R I  NT , THIN 

7  FORMAT  (10F8.3) 

IF    ( NKOLS    . LE.    0)    GO   TO  8 
READ    (5,6)    (KOL(I), 1=1, NKOLS) 

8  READ   (5,7)    AGMRCH , BFMRCH , B FS AL V , C OMCU , COMBF , BF PCT , C FPCT , GNTR 
READ    (5,7)    BDPRI, CFPRI 

CRE4TE    A   SERIES   OF   CONSTANT   OR   VARIABLE  PRICES 

IF    (CFPRI    .NE.    0.0)    GO   TO  9 
READ    (5,7)    (PRICFt I), 1=1,150) 
GO   TO  11 

9  DO    10  1=1,150 

10  PRICFt I )    =  CFPRI 

11  IF  (BDPRI  .NE.  0.0)  GO  TO  12 
READ  (5,7)  (PRIBDI I ) , 1=1, 150) 
GO   TO  14 

12  DO    13  1=1,150 

13  PRIBDI I )    =  BDPRI 

14  RETURN 
END 

SUBROUTINE  YIELD 

COMMON  ACC ST, AGEO, AGMRCH, AN BDF( 181),ANCUV(181),ANNET,BATCH(3),BDFC 
1 ( 180) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 180) , C STAC , C ST VL , CUCST , CUT AGE 
2,DBH0,DEN0,DESCR(5 ) , DLE V , FMRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V AL B , GS V ALC , G VL B F , 
3GVLCU, I  ACRE! 180) , I ALCUT , I  GAME , ISUMI 18) , I TE ST, I VAR ( 26 , 1 50 ) , I  YEAR, 
4K0L(6) ,LAND,LAST,MALCUT( 10) .MOLD, NACRE ( 1 80 ), NG AME , NKOLS , Nl , NONS T K , 
5N0YRS, PRET, PRIBDI 150) ,PRICF(150),PRIDIV( 1 0 ) , RE TRN , ROTA , RATE , GROW , 
6SITE,SUMM(6,25, 10 ) , TCOST , SHELT , YRLOS , F I NL , C LOS S , CTHN , C PLT , CFPC T , 
7VAR(14,150),VLBF,VLCU,ACAGE(1000) , C YCL , KOUNT , DEFOR , ANUL , MIX,  I  PL  NT, 
8RI NT, THIN, BFMRCH, BFSALV, COMCU, COMBF, BFPCT,GNTR 

DIMENSION   TABLK  81 ) ,TABL2<4,18) ,TABL3(192) , TABL4 (156) 

BDFT   =  0.0 

CFMT   =  0.0 

JBDFC   =  0 

JBDFO  =  0 

JBDFT   =  0 

JCFMC   =  0 

JCFMO   =  0 

JCFMT   =  0 

DO   1  1=1,180 

BDFC ( I )  =0.0 

BDFO ( I )    =  0.0 

CFMC( I )    =  0.0 

1  CFMO( I )    =  0.0 

READ   HEIGHTS, DENSITIES, VOLUME   CONVERSION  FACTORS 

READ    (5,2)    ( TABL1 IK) ,K=1 ,81 ) 

2  FORMAT  (21F3.1) 

READ    (5,3)    ( (TABL2IK,L),K=1,4) ,L=1,18) 

3  FORMAT  (25F3.1) 

READ    (5,4)    1TABL3(K),K=1, 192) 

4  FORMAT  (24F3.3) 

READ    (5,5)    (TABL4IK) ,K=1 , 156) 

5  FORMAT  (26F3.2) 
Nl   =  AGEO 

N  =  AGEO 

BASO   =   DENO   *   0.0054542    *   DBHO   *  DBHO 
OBTAIN  HTSO 
TEN   =  10.0 

ISITE   =    (SITE/TEN  -   3.0)    +  0.01 
I AGEO  =    (AGEO/TEN)    +  0.01 
HTSO   =   TABL2I ISITE, IAGEO) 

TOTO   =    (0.4047   *   BASO   »   HTSO)    +    (25.5970   »   DBHO)    -  191.6433 

COMPUTE   MERCH.    CU.FT.    IF   DBH    IS    AT   LEAST   5.0  INCHES 

IF    (DBHO    .LT.    5.0)    GO   TO  6 

I DBHO  =    ((DBHO  -   5.0)    *   10.0)    +  1.01 

XDBHO   =   TABL3I I DBHO ) 

CFMO(N)    =   TOTO   •  XDBHO 

COMPUTE   BD.FT.    IF   DBH   IS    AT   LEAST   8.0  INCHES 

IF    (DBHO    .LT.    8.0)    GO   TO  6 

JOBHO   =    ((DBHO  -   8.0)    *    10.0)    +  1.01 

YDBHO  =   T  ABL4 ( J  DBHO ) 

BDFO ( N )    =   TOTO   *  YDBHO 

COMPUTE   DBH   AFTER    INITIAL  THINNING 

6  DO   11  J=l,100 

PDBHE   =   0.95462*AL0G10(DBH0)    -   0 . 10640* ALOG 1 0 ( PRE T )    +  0.26959 

DBH  E   =    10.0   »*  PDBHE 

IDBHE   =   DBH  E   *    10.0   +  0.5 

DBHE   =  IDBHE 

DBHE   =  DBHE/TEN 

DENE   =   DENO   »  (PRET/100.0) 

BASE   =    (0.0054542    »   DBHE    »   DBHE)    *  DENE 

BREAK   =   49.9   »  THIN/80.0 

IF    (BASE    .  GT.    BREAK)    GO   TO  7 

D8HP   =    (80.0/THIN)    *    (0.08733  »   BASE)    +  0.92247 
GO  TO  8 

7  DBH  P   =    (80.0/THIN)    *  (0.10938 

8  IDBHP   =   DBHP   •    10.0   +  0.5 
DBHP   =  IDBHP 
DBHP   =  DBHP/TEN 
IF    (DBHP  -  DBHE)  9,12,10 

9  PRET   =   PRET   +   PRET   *  0.02 
GO   TO  11 


BASE)    -  0.17858 


10  PRET   =   PRET   -  PRET   »  0.02 

11  CONTINUE 

12  DBHT   =  DBHE 

COMPUTE   VALUES    AFTER   INITIAL  THINNING 

JDBHT   =    ((DBHT  -   2.0)    »    10.0)    +  1.01 

SQFT   =   TABL1 ( JDBHT ) 

BAST   =    (THIN/80.0)    *  SQFT 

ENTER   LOOP   FOR   ALL   REMAINING  COMPUTATIONS   AND  PRINTOUT 


DO   26  1=1,100 

DENT   =   BAST/IO. 0054542 

HTST   =  HTSO 

TOTT   =    (0.4047   »   BAST  • 


f  DBHT 
HTST) 


DBHT  ) 

( 25.5970 


COMPUTE   MERCH.    CU.FT.    IF   DBH   IS   AT   LEAST   5.0  INCHES 

IF    (DBHT    .LT.    5.0)    GO  TO  13 

I DBHT  =    ((DBHT   -   5.0)    »    10.0)    +  1.01 

XDBHT   =   TABL3 ( I DBHT ) 

CFMT   =   TOTT   *  XDBHT 

COMPUTE   BD.FT.    IF   DBH    IS    AT   LEAST   8.0  INCHES 

IF    (DBHT    .LT.    8.0)    GO   TO  13 

JDBHT   =    ((DBHT   -  8.0)    *   10.0)    +  1.01 

YDBHT   =   T  ABL4 ( JDBHT ) 

BDFT   =   TOTT   »  YDBHT 

CHANGE   MODE   AND   ROUND  OFF   FOR  PRINTING 


13  JDENO 

jhtso 

JTOTO 
JBASO 
JCFMO 
JBDFO 
JBDFO 
JDENT 
JHTST 
JTOTT 
JDENC 
JCFMT 
CFMT 


DENO 
HTSO 

TOTO  +  0.5 
BASO   +  0.5 
CFMO(N)    +  0.5 
( BDFO ( N ) /TEN )    +  0.5 
JBDFO  *  10 
DENT 
HTST 

TOTT   +  0.5 
JDENO  -  JDENT 
CFMT   +  0.5 
JCFMT 

IF    (JCFMT    .GT.    JCFMO)    JCFMO  =  JCFMT 

CFMO(N)    =  JCFMO 

JBDFT  =    (BDFT/TEN)    +  0.5 

JBDFT   =   JBDFT   «  10 

BDFT   =  JBDFT 

BDFT  =   BDFT   •  .001 

IF    (JBDFT    .GT.    JBDFO)    JBDFO  =  JBDFT 

BDFO(N)    =  JBDFO 

BDFO(N)    =   BDFO ( N )    *  .001 

JBAST  =   BAST   +  0.5 

JBASC   =   JBASO  -  JBAST 

JTOTC   =   JTOTO  -  JTOTT 

JCFMC   =   JCFMO  -  JCFMT 

IF    (JCFMC    .LE.    0)    JCFMC   =  0 

CFMC(N)    =  JCFMC 

JBDFC   =   JBDFO  -  JBDFT 

IF    (JBDFC    .LE.    0)    JBDFC   =  0 

BOFC(N)    =  JBDFC 

BDFC(N)    =   BDFC(N)    »  .001 

IF    (I    .GE.    2)    GO   TO  19 

WRITE   HEADINGS    FOR   YIELD  TABLE 

WRITE    (6,14)  SITE,CYCL,DLEV 

14  FORMAT  ( 1H1,////,28X,81HYIELDS  PER  ACRE  OF  MANAGED,  EVEN-AGED  STAN 
IDS  OF  PONDEROSA  PINE  IN  THE  BLACK  HILLS/1H  ,35X,10HSITE  INDEX, F3.0 
2,1H, , F4.0, 19H-YEAR   CUTTING  C YC L E , 1 H , , 1 5H   DENSITY   LEVEL  ,F4.0,///) 

WRITE  (6,15) 

15  FORMAT  (1H0,25X,38HENTIRE  STAND  BEFORE  AND  AFTER  THINNING, 28X,26HP 
1ERI0DIC   CUT   AND  MORTALITY) 

WRITE  (6,16) 

16  FORMAT  ( 1 HO, 9X.5HST AN D,10X,5HBASAL,3X,7H AVER AG E,2X,7H AVERAGE, 3X.5H 
1T0TAL,3X,9HMERCHANT-,3X,9HSAWTIMBER,9X,5HBASAL,4X,5HT0TAL,3X,9HMER 
2CHANT-,3X,9HSAWTIMBER) 

WRITE  (6,17) 

17  FORMAT  (1H  ,10X,3HAGE,4X,5HTREES,3X,4HAREA,4X,6HD.B.H.,3X,6HHEIGHT 
1, 2 X,6H VOLUME, 2X, 11H ABLE  VO LUME , 4X , 6H VOL UME , 3X , 5HTREE S , 3 X , 4H ARE  A , 3 X 
2,6HV0LUME,2X, 11HABLE   VOLUME , 4X , 6H VOL UME ) 

WRITE  (6,18) 

18  FORMAT  (1H  ,8X,7H( YEARS) , 3X , 3HN0 . , 3 X , 6HSQ . F T . , 4X , 3H I N . , 6X , 3HFT . , 4X 
1,6HCU.FT.,5X,6HCU.FT. , 8X , 3HMBF , 5X , 3HN0 . , 3X , 6HSQ . FT . , 2X,6HCU.FT.,5X 
2.6HCU. FT. ,8X,3HMBF ) 

19  WRITE    (6,20)    AGEO , JDENO , J B A  SO , DBHO , JHTSO , JTOTO , C FMO ( N ), BDFO ( N ) 

20  FORMAT  (1HO,9X,F4.0,4X,I5,2X,I4,5X,F5.1,5X,I3,4X,I5,6X,F5.0,6X,F6. 
13) 

IF    (AGEO. GE. ROTA)    GO  TO  27 

WRITE  (6,21)  AGEO, JDENT, JBAST, DBHT, JHTST, J TOTT, CFMT, BDFT, JDENC, JBA 
ISC, JTOTC, CFMC(N) ,BOFC(N) 

21  FORMAT  ( 1H  , 9X , F4 . 0 , 4X , I  5 , 2X , 1 4, 5X , F 5 . 1 , 5X , I  3 , 4X , I  5 , 6X , F5 . 0 , 6X , F6 . 
13,4X,I5,3X,I3,5X,I4,6X,F4.0,8X,F5.3) 

COMPUTE   VALUES   FOR  EACH   PERIOD.    THIN   AS  SPECIFIED 


IK   =  CYCL/RINT 
DO   23  L-ltlK 
AGEO  =   AGEO  +  RINT 
N  =  AGEO 

IF    (AGEO   .    GT    .    ROTA)    GO   TO  27 
DBHO  =   1.0097»DBHT   +  0.0096*SITE 
I DBHO  =   DBHO   ♦   10.0  +  0.5 
DBHO  =    I DBHO 
DBHO   =  DBHO/TEN 


1. 5766* ALOGlOt BAST) 


-   17  - 


ADO  RANDOM  ELEMENT  TO  PREDICTED  DBHOtIF  DESIRED 

IF    (GNTR    .GT.    1024.0)  GO   TO  110 

IDIV   =    (17.0   »   GNTR   +  3.01/1024.0 

NGNTR  =  GNTR 

GNTR   =    (17   ♦   NGNTR   +    3)    -    1024  *  IDIV 

IF    (GNTR    .GT.    1000.0)  GO   TO  100 

IF    (GNTR   . L T .   0.0)    GO  TO  100 
Al    =  GNTR/100.0 


Al  -   0.0523   *   A2  -  0.0063 


IRES  =  RES 
IRES  =  RES 


0.5 
0.5 


DBHO 
0.01 


A2   =   Al    *  Al 
RES   =  0.9565 
1*A2   -  3.3009 
IRES   =  RES 
IF    (RES    . LT .  0.0) 
IF    (RES    .GT.  0.0) 
AOJ    =  IRES 

DBHO   =   DBHO   +  ADJ/TEN 
110   DENO   =  DENT 

BASO  =  DE  NO  *  (0.0054542 
ISITE  =  (SITE/TEN  -  3.0) 
IAGEO  =  (AGEO/TEN)  +  0.01 
HTSO  =  TABL2I ISITE, IAGEO) 
TOTO  =  0.4047  «   BASO  *  HTSO 


COMPUTE   MERCH.    CU.FT.    IF   OBH    IS   AT   LEAST   5.0  INCHES 

IF    (DBHO   .LT.    5.0)    GO  TO  22 

IDBHO  =    ((DBHO  -   5.0)    »    10.0)    +  1.01 

XDBHO   =   TABL3 ( IDBHO) 

CFMO(N)    =   TOTO   *  XDBHO 

COMPUTE   BD.FT.    IF   DBH    IS    AT   LEAST   8.0  INCHES 

IF    (DBHO    .LT.    8.0)    GO   TO  22 

JDBHO   =    ( (DBHO  -   8.0)    *    10.0)    +  1.01 

YDBHO   =   T ABL4 ( JDBHO ) 

BDFOIN)    =   TOTO   *  YDBHO 

IF    ( L    .EQ.    IK)    GO   TO  24 

WRITE   VALUES    FOR   END   OF   PERIOD    IF   THINNING   NOT  DUE 


KDENO  =  DENO 
KHTSO  =  HTSO 
KBASO  = 
KTOTO  = 
JCFMO  = 
CFMO(N) 
JBDFO  = 
JBDFO  = 
BDFOIN) 


1.5 


0.5 


;toto,cfmo(n) .bdfoin) 


BASO  +  0.5 
TOTO  +  0.5 
CFMO(N)    +  C 
=  JCFMO 
(BDFO(N)/TEN) 
JBDFO  *  10 
=  JBDFO 
BDFO(N)    =   BDFOIN)    »  .001 

WRITE    (6,20)    AGEO, KDENO, KBASO, DBHO, KHTSO, 
OBHT   =  DBHO 
BAST   =  BASO 
23  CONTINUE 


INCREASE   DBH   AS   RESULT   OF   THINNING  AND   COMPUTE   PO ST- TH I NN I NG 
VALUES 

24  DBHT   =    DBHO   +  0.4 

IF    (DBHT    . GE .    10.0)   GO  TO  25 

JDBHT   =    ((DBHT   -   2.0)    •    10.0)    +  1.01 

SOFT   =   T ABL 1 ( JDBHT ) 

BAST   =    (DLEV/80.0)    «  SOFT 

GO   TO  26 

25  BAST  =  DLEV 

26  CONTINUE 

27  RETURN 
END 

SUBROUTINE  ANVOL 

COMMON   ACCST.AGEO.AGMRCH, ANBDF( 1811, ANCUV ( 181 ) , ANNE T , BATCH ( 3 ), BDFC 
1(18  0)  ,BDFO(  180),BFCST,CFMC(180),CFM0(180) , C STAC , C STVL , CUC ST , CUT AGE 
2,DBH0,DEN0,DESCR(5)  ,  DL E V , F MRCHD ( 10) ,GMNAM( 3) , GS VALB , GSVALC , GVLBF , 
3GVLCU, I  ACRE! 180) , IALCUT, I  GAME, I  SUM! 18),ITEST,IVAR(26,150),IYEAR, 
4K0LI6) , LAND, LAST, MALCUTI 10) , MOLD, NACRE (180) , NG AME , NKOL S , Nl , NONSTK , 
5N0YRS.PRET, PRIBDt 150),PRICF(150),PRI0IV( 10) ,RETRN, ROTA, RATE, GROW, 
6SITE,SUMM(6,25, 10 ) , TCOST , SHELT , YRLOS , F I NL , CLOSS , C THN , CPLT , CFPCT , 
7VAR(14,150) ,VLBF,VLCU,ACAGE( 1000)  , C YCL , KOUNT , 0 E FOR , ANUL , M I  X , IPLNT, 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT, GNTR 

IROT   =  ROTA 

INT   =   R I  NT 

NVUL   =    (  (  IROT   -   N  1  )  /  I  NT  )    +  1 
K   =   NVOL   -  1 


WRITE   CUBIC   FEET   PER   ACRE   FOR   EACH  YEAR 

WRITE    (6,5)    K, ( ANCUV(NN) ,NN=1, 10) 

5  FORMAT    (1H    , I  20 , F 1 3. 1 , 9F 1 0 . 1 ) 
IJ   =    IR0T/10  -  1 

DO   6   J  =  l , I  J 

NN  =    10    *   J   +  1 

WRITE  (6,5)  J, ANCUV (NN) , ANCUV (NN+1) , ANCUV (NN+2 ) , ANCUV (NN+3) , ANCUV ( 
1NN+4) ,ANCUV(NN+5) , ANCUV ( NN+6 ) , ANCUV ( NN+7 ) , ANCUV ( NN+8 ) , ANCUV ( NN+9 ) 

6  CONTINUE 

J  =  IROT/10 

ANCUVlIROT+l)    =  CFMO(IROT) 
WRITE    (6,5)    J , ANCUVI IROT+1 ) 

WRITE   BD.FT.    PER  ACRE  FOR   EACH  YEAR 

WRITE  (6,7) 

7  FORMAT    ( 1H0,///,55X,23HTH0USANDS   OF   BOARD  FEET,//) 
WRITE    (6,8)    K, ( ANBDF I NN) ,NN=1, 10) 

8  FORMAT    I  1H   , I  20 , F 1 3 . 3 , 9F 1 0 . 3 ) 
DO   9   J= 1  ,  I  J 


NN 


10 


WRITE  (6,8)  J, ANBDF (NN) , ANBDF (NN+1) , ANBDF (NN+2) , ANBOF ( NN+3 ), ANBDF ( 
1NN+4) , ANBDF ( NN+5 ) , ANBDFI NN+6), ANBDF (NN+7), ANBDF(NN+8),ANBDF( NN+9) 

9  CONTINUE 
J    =  IROT/10 

ANBDFI IROT+1)    =   BDFO ( IROT) 
WRITE    (6,8)    J , ANBDF ( IROT+1 ) 

PROVIDE   FOR   ANY  ACRES  BEYOND  CUTAGE   LEFT  UNTHINNED  FOR  A  FEW  YEARS 

JCYCL   =  CYCL 

IMI    =   JCYCL  -  1 

DO   10    I=N1, IROT, JCYCL 

DO    10  J=1,IMI 

NX  =   I   +  J 

IF    (NX    .GE.    180)    GO   TO  11 
BDFC ( NX )    =   BOFCI I) 

10  CFMC(NX)    =   CFMCI  I  ) 

11  RETURN 

END 

SUBROUTINE  INPUT2 

COMMON   ACCST.AGEO.AGMRCH, ANBDFI 181) , ANCUVI 181) , ANNE T , BATCH  (  3 ),  BDFC 
1  ( 180) . BDFO ( 180) , 8FCST, CFMCI 180) ,CFMO( 180 ) , C ST AC , C ST VL , C UC ST , CUT AGE 
2,DBH0,DEN0,DESCR( 5) , DLE V , F MRCHD I  1 0 ) , GMNAM I  3 ) , GS V A  LB , GSVALC , GVLBF , 
3GVLCU,IACRE(180) , IALCUT, IGAME,ISUM(18),ITEST,IVAR(26,150),IYEAR, 
4K0L( 6) , LAND, LAST, MALCUTf 10 ), MOLD, NACRE ( 180 ),NG AME, NKOLS.Nl, NONSTK, 
5N0YRS,PRET,PRIBD( 150 ) , PR  I CF ( 1 5 0 ) , PR  I D I  V ( 10 ) , RE T RN , ROT A , RATE , GROW , 
6SITE,SUMM(6, 25, 10), TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT, CFPCT, 
7VAR( 14, 150)  ,VLBF,VLCU, AC AGE ( 1000) , C YCL , KOUNT , D EF OR , ANUL , M I  X , I  PL  NT , 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT,GNTR 

SET    INITIAL   VALUES   OF  ZERO 

CSTAC  =  0.0 
CSTVL   =  0.0 
IALCUT   =  0 
KOUNT   =  1 
LAST   =  0 
RETRN   =  0.0 
YRLOS   =  0.0 
DO    1  1=1,180 
I  AC  RE (  I  )    =  0 

1  NACRE ( I )    =  0 
DO   2  1=1,10 
FMRCHD ( I )    =  0.0 
MALCUT(I)    =  0 

2  PRIDIVI I )    =  0.0 
DO   3  1=1,26 

DO   3  J=l,150 

3  IVARI I , J)    =  0 
DO   4    1=1, 14 
DO  4  J=l,150 

4  VARI I, J)    =  0.0 

READ   VALUES   THAT   DO   NOT   CHANGE   WITHIN   A  GAME 

READ    (5,5)    ( GMNAM ( I ) , I = 1 , 3 ) 

5  FORMAT  (3A8) 

READ    (5,6)    LAND, MOLD, NONSTK, KAREA, IPLNT 

6  FORMAT  (2014) 

IF    (KAREA   . EO.    0)   GO  TO  8 


INTERPOLATE    BETWEEN   VALUES   FROM   YIELD  TABLE 


READ    IN   EQUAL   AREAS    IN   ALL   AGE  CLASSES 


(CFMO(N+INT)-CFMO(N)+CFMC(N) ) 
(BDFO(N+INT)-BDF0(N)+BDFC(N) ) 


DO    1  L=1,K 
DO    1  J=1,INT 

NN   =   J    +   Nl   +    (L   -   1)    *  INT 
RJ   =   J   -  1 

N   =   Nl   +    (L  -   1)    »  INT 

ANCUV (NN)  =CFMO(N)-CFMC(N)+(RJ/RINT)i 
ANBDF (NN)    = BDFO ( N ) -BDFC ( N )  +  ( R J / R I  NT ) < 
1  CONTINUE 

WRITE    TABLE  HEADINGS 


WRITE    (6,2)    SITE, CYCL, THIN, DLEV 

2  FORMAT  ( 1H1,////,41X,51HGR0WING  STOCK  OF  MANAGED  BLACK  HILLS  PONDE 
1R0SA  PINE/1H  ,47X,10HSITE  I NDE X , F 3 . 0 , 1H , , F4. 0 , 1 9H- YE AR  CUTTING  CYC 
2LE/1H    ,53X, 14H DENSITY   LE VEL- , F4 . 0 , 1 X , 3H AND, F4 . 0 ) 

WRITE  (6,3) 

3  FORMAT    ( 1HO,43X,44HVOLUMES   PRESENT   PER   ACRE   AT    END  OF   EACH  YEAR) 
WRITE  (6,4) 

4  FORMAT  ( 1H  , 54X,23H MERCHANTABLE  CUBIC  F EET/ 1H0 , 64X , 4HYE AR/ 1 H  ,14X, 
16HDECADE,9X,1H0,9X, 1H1.9X, 1H2 , 9X , 1 H3 , 9X , 1H4 , 9X , 1 H5 , 9X , 1H6 , 9X , 1H7 , 9 
2X, 1H8.9X, 1H9,// ) 


NDX   =   MOLD  +  1 

DO   7    1=1, NDX 

I  ACRE ( I  )    =  KAREA 

ADJUST  NUMBER  OF  ACRES  IN  OLDEST  CLASS  IF  TOTAL  AREA  NOT  MULTIPLE 
OF  KAREA 

KDIFF   =    LAND  -   NDX   *   KAREA  -  NONSTK 
I  ACRE ( NDX )    =   I  ACRE ( NDX )    +  KDIFF 
IACRE(l)    =   IACRE(l)   +  NONSTK 
GO  TO  10 

READ   IN  UNEQUAL   AREAS    IN  AGE  CLASSES 


READ   LIMITATIONS   ON  CUT 


10  READ    (5,11)    ( PRIDIVI I ), 1=] 

11  FORMAT  (10F8.3) 
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READ    (5,9)    (MALCUTI I ), 1=1,10) 
READ    (5,11)    (FMRCHDI I ), 1=1,10) 

READ  15,11)  SHELT, RATE, CPLT, CTHN, CLOSS, ACCST, CUCST, BFCST, GROW, FINL 
READ    (5,12)    DEFOR , ANUL 

12  FORMAT  (2F8.5) 
DEF0R1   =   DEFOR   »  100.0 

PRINT   CONDITIONS   OF  SIMULATIONS 

WRITE  (6,13) 

13  FORMAT    ( 1H1 , // , 46X.26HALTERNAT I VES   FOR   THIS  GAME) 
WRITE    (6,14)    (BATCH! I ), 1=1,3) 

14  FORMAT    (1H    ,45X,7HBATCH  ,3A8) 
WRITF    (6,15)  ITEST 

15  FORMAT    (1H    , 45X , 4HTE ST , I  4 ) 
WRITE    (6,16)    ( GMNAM ( I ) , 1=1 ,3 ) 

16  FORMAT    ( 1H   ,45X,6HGAME  ,3A8) 
WRITE    (6,17)    (DESCRI I ) , 1=1 ,5) 

17  FORMAT    (1H   , 45X , 5 A8 , / / / / ) 
WRITE    (6,18)  NOYRS 

18  FORMAT    (1H    , 45 X , 24HNUMBER  OF   YEARS   PER   GAME , 14,////) 
WRITE    (6,19)    ( PRIDIVI I ) , 1=1, 10) 

19  FORMAT    (1H    .15HCRITICAL   PR  I C E S , 1 2 X , 1 0F9 . 2 ) 
WRITE    16,20)    ( MALCUT (I), 1=1, 10) 

20  FORMAT    (1H    , 13H ALLOWABLE   C UT , 1 1 X , 1 0  I  9 ) 
WRITE    (6,21)    ( FMRCHDI I ), 1=1 , 10) 

21  FORMAT    ( 1H    , 19HM I N I  MUM   CUTTING   AGE , 5 X , 1  OF  9 . 0 , / / / / ) 
WRITE    (6,22)  LANO 

22  FORMAT    (1H    ,23HACRES    IN  WORKING   CIRCLE, 13X, I4,25X,27HC0STS   IN  FIRS 
IT   YEAR   OF  GAME) 

WRITE    (6,23)  ACCST 

23  FORMAT    (1H    ,69X,17HPER   ACRE    (ANNUAL), 8X,F9. 2) 
WRITE    (6,24)  CUCST 

24  FORMAT    (1H    .3RHMINIMUM   VALUES   FOR    INCLUSION    IN   TOTALS, 31X,25HPER  1 
100   CU.    FT.    HARVESTED, F9. 2) 

WRITE    (6,25)  AGMRCH.BFCST 

25  FORMAT    ( 1H    ,4X,22HAGE,    FOR   GROWING   STOCK , 1 1 X , F 3 . 0 , 29X , 13HPER  M  BD. 
1   FT., 12X.F9.2) 

WRITE    (6,26)    BFMRCH, CTHN 

26  FORMAT    (1H    ,4X,28HM   BD.    FT.,    FOR   GROWING   STOCK , 5X , F5 . 1 , 27X , 1 3HTH I N 
1   ONE  ACRE.12X.F9.2) 

WRITE    (6,27)  COMCU.CPLT 

27  FORMAT    ( 1H    ,4X,27HCU.    FT.,    FOR   COMMERCIAL   CUT , 5 X , F4. 0 , 29X , 14HPL A  NT 
1   ONE  ACRE.11X.F9.2) 

WRITE    (6,28)  COMBF.CLOSS 

28  FORMAT    (1H    ,4X,29HM   BD.    FT.,    FOR   COMMERCIAL   CUT , 4X, F 5 . 1 , 27X , 1 9HC LE 
1ANUP   OF   ONE  ACRE.6X.F9.2) 

WRITE    (6,29)  BFSALV.RATE 

29  FURMAT    ( 1H    ,4X,22HM   BD.    FT.,    FOR   S AL V AGE , 1 1 X , F 5 . 1 , 23X , 2 5HR ATE   OF  I 
INCREASE    IN  COSTS , 4X , F9 . 2 , / / ) 

WRITF    (6,30)  IPLNT 

30  FORMAT    (1H   , 22H ACRES   PLANTED  ANNUALLY, 14X, 14, 25X,35HRELATIVE  VALUE 
1   OF    INTERMEDIATE  CUTS) 

WRITE    (6,31)  DEF0R1,CFPCT 

31  FORMAT    (1H    , 30HPERCENT   OF   ACRES   LOST   ANNU ALL Y , 6X , F8 . 3 , 2 5X , 23HSTUMP 
1  AGE   PRICE,    CU.  FT.,2X,F9.2) 

WRITE    (6,32)    SHELT, BFPCT 

32  FORMAT    ( 1H    ,24HM  BD.    FT.    IN   SHELT ERWOOD , 1 3X , F5 . 1 , 27X , 2 3HSTUMP AGE  P 
lRICfc,    BO.  FT.,2X,F9.2) 

WRITE    (6,33)  FINL 

33  FORMAT    ( 1H    , 1 9HREGENER AT  I  ON   P ER I OD, 1 8X , F 5 . 1 , // ) 
WRITE    (6,34)  ANUL 

34  FORMAT    ( 1H    , 29HPSEUD0RAND0M  NUMBER  GENERATOR, 5X,F8.1) 
WRITE    (6,35)  GNTR 

35  FORMAT    (1H  ,34X,F8.1) 
RETURN 

END 

SUBROUTINE  AREAS 

COMMON   ACCST, AGEO, AGMRCH, ANBDF ( 181) , ANCUVI 181 ) , ANNE T , BATCH ( 3 ) , BDFC 
1(180) ,BDF0(180) ,BFCST,CFMC( 180) ,CFMO(  180) , C ST  AC , C ST VL , C UC ST , CUTAGE 
2 , DBHO , DENO, DESCR(5),DLEV, FMRCHDI 10) , GMNAM ( 3) , GS VA LB , GS V ALC , G VLBF , 
3GVLCU, I  ACRE (180) , I ALCUT, I  GAME, ISUMt 18) , I TES T , I  V AR ( 26 ,  150) , I  YEAR, 
4K0L (6) , LAND, LAST, MALCUT ( 10 ) , MOLD, NACRE ( 180 ) , NGAHE , NKOLS , Nl , NONSTK, 
5N0YRS, PRET.PRIBDt 150),PRICF(150),PRIDIV(10) , RE TRN , ROT A , RATE , GROW , 
6SITE,SUMM(6,25,10) , TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT,CFPCT, 
7VAR(  14, 150) ,VLBF,VLCU, AC  AGE ( 1000), CYCL,KOUNT, DEFOR, A NUL, MIX, IPLNT, 
8RI NT, TH I N, BFMRCH, BFS AL V , COMCU , COMBF , BFPCT,  GNTR 

GSVALB   =  0.0 

GSVALC   =  0.0 

GVLBF   =  0.0 

GVLCU  =  0.0 

00    1  1=1,18 

1  ISUM( I )    =  0 
DO   2  1=1,1000 

2  ACAGE ( I )    =  0.0 

CONVERT    ACRES    IN   EACH    IACRE(I)    TO    INDIVIDUAL  ACRES 

JK  =  0 

DO   30  J=l,180 

IF    (JK    .GE.    LAND)    GO   TO  4 

IF    ( IACREI J)    .LE.    0)    GO   TO  30 


KL 


JK 


JK   =   JK   +  IACREIJ) 
DO   3  I=KL,JK 
NAC   =   LAND  +1-1 
ACAGE I NAC )    =   J   -  1 
3  CONTINUE 
30  CONTINUE 


DO   7   K=l , LAND 

IF    ( ACAGE ( K )    .LE.    179.0)    GO   TO  6 
WRITE  (6,5) 

FURMAT  ( 1H1 ,47X,38HY0U  WENT  BEYOND  AGE  LIMIT  OF  179  YEARS) 
IYEAR   =   NOYRS   -  1 


RETURN 

6  LM   =   ACAGE ( K )    +  1.0 

7  NACRE ( LM )    =   NACRE ( LM )    +  1 

COMPUTE   TOTAL   ACREAGE   BY    10-YEAR   AGE  CLASSES 


8  ISUM(I)    =    ISUM(I)    +  NACRE ( NS ) 

COMPUTE   GROWING  STOCK  VOLUME. 

USE   CU.FT.    IF   VOLUME    IS   LESS   THAN  BFMRCH 

DO    10   M= 1 , L  AND 

IF    ( ACAGE ( M )    .LT.    AGMRCH)    GO  TO  10 
I  AG   =    AC  AGE ( M )    +  1.0 

IF    ( ANBDF ( I  AG )    .GE.    BFMRCH)    GO   TO  9 
GVLCU  =   GVLCU   +   ANCUVI  I  AG ) 
GO   TO  10 

9  GVLBF   =   GVLBF   +  ANBDFIIAG) 
10  CONTINUE 


COMPUTE  IN 

IVAR(7,1) 
IVAR(8,1) 
IVARI 9, 1 ) 
IVARUO.l) 
IVARI 11,1) 
VAR I  1 , 1 )  = 
VARI2.1)  = 
GSVALB  = 
GSVALC  = 
VARI 13, 1 
VARI 14, 1 
DO    11  1= 


I  AL   NON-ZERO  VALUES   FOR  0UTPUT2 

GVLCU   +  0.5 
GVLBF   +  0.5 
IVAR(5,1)    +  IVAR(7,1) 
I VAR (6,1)    +    IVARI 8,1) 
)    =  NONSTK 
=   PRICFI 1) 
=   PRIBOI 1 ) 

GVLBF   »    (PRIBD(l)    -  BFCST) 
(GVLCU/100.0 )    *    (PRICF(l)    -  CUCST) 
=   GSVALC   +  GSVALB 
=   VARI 13,  1  )    +   VARI 12,  1  ) 
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WRITE   HEADINGS   FOR   TABLE    OF    INITIAL   DISTRIBUTION   OF    ACRES   BY  AGE 
WRITE  (6,12) 

12  FORMAT    ( 1H1,////,38X,36HINITIAL   DISTRIBUTION  OF   ACRES   BY  AGE) 
WRITE    (6,13)    (BATCH!  I  ), 1  =  1,3) 

13  FORMAT    ( 1H    ,45X,7HBATCH  ,3A8) 
WRITE    (6,14)  ITEST 

14  FORMAT    (1H    , 45X , 4HTE S T ,  I  4 ) 
WRITE    (6,15)    (GMNAM (I), 1=1, 3) 

15  FORMAT    (1H    ,45X,6HGAME  ,3A8) 
WRITE    (6,16)    (DESCRI I), 1=1, 5) 

16  FORMAT    ( 1H    ,45X,5A8 ) 
WRITE  (6,17) 

17  FORMAT    (1H    ,45X,16HYEAR   WITHIN   GAME , 3X , 1H0 , /// J 
WRITE  (6,18) 

18  FORMAT    ( 1H    , 55X , 9HAGE ( YE AR ) ) 
WRITE  (6,19) 

19  FORMAT    (1H    ,4X, 11HAGEI DECADE) ,7X, 1H0,7X,1H1,7X,1H2,7X,1H3,7X,1H4,7 
1X,1H5,7X,1H6,7X,1H7,7X,1H8,7X,1H9,10X,5HT0TAL,//) 


WRITE  NUMBER  OF  ACRES  IN  EACH 
10-YEAR  CLASSES 


-YEAR   AGE   CLASS   AND   TOTALS  OF 


10 


IK 


WRITE  (6,20)  IK, NACRE ( NN ), NACRE (NN+1 ), NACRE (NN+2) , NACRE (NN+3),NACR 
1EINN+4) , NACRE (NN+5) , NACRE (NN+6) ,NACRE(NN+7),NACRE(NN+8),NACRE(NN+9 
2 ) , I  SUM ( J ) 

20  FORMAT    ( 1H    ,111, 5X, 1018, 115,/) 

21  CONTINUE 
RETURN 
END 

SUBROUTINE  YEARS 

COMMON   ACCST, AGEO, AGMRCH, ANBDF (181), ANCUVI 181 ),ANNET, BATCH (3), BDFC 
1(180), BDF0I180) , BFCST, CFMC (180) ,CFMO( 180) ,CST ACCST VL, CUCST, CUT AGE 
2,DBH0,DEN0,DESCR(5)  , DL E V , FMRCHD I  1 0 ) , GMN AM  I  3 ) , GS VA L B , GS VALC , G VL BF , 
3GVLCU, I ACRE( 180) ,  I ALCUT,  I  GAME , I  SUM ( 18) , I  TEST , I VAR ( 26 , 1 50 ) , IYEAR, 
4K0L ( 6), LAND, LAST, MALCUT (10), MOLD, NACRE ( 180 ),NGAME, NKOLS, Nl, NONSTK, 
5N0YRS,PRET,PRIBD(  1 50 ) , PR  I CF ( 150 ) , PR  I D I V ( 1 0 ) , RE TRN , ROT A , RA TE , GROW , 
6SITE,SUMM(6, 25,10), TCOST, SHELT, YR LOS, FINL, CLOSS, CTHN, CPLT,CFPCT, 
7VAR(14,150),VLBF,VLCU,ACAGE(1000) , C YCL , KOUNT , DEFOR , ANUL , M I  X , I  PL  NT , 
8R I  NT, THIN, BFMRCH, BFS ALV, COMCU, COMBF, BFPCT, GNTR 

GVLBF   =  0.0 

GVLCU   =  0.0 

LOSS   =  0 

NPLNT  =  0 

RETHV   =  0.0 

RETTH   =  0.0 

SCLOSS   =  0.0 

SCPLT   =  0.0 

SCTHN  =  0.0 

VLBF   =  0.0 

VLCU   =  0.0 

JCYCL   =  CYCL 

IYRM  =    IYEAR   +  1 


MAKE 


SCHEDULED   ANNUAL  PLANTING 


IF    (NONSTK    .  EQ .    0)    GO  TO  1 
NPLNT   =  IPLNT 

IF    (NPLNT    .GT.    NONSTK)    NPLNT   =  NONSTK 

NONSTK   =  NONSTK  -NPLNT 

APLT   =  NPLNT 

SCPLT   =   APLT   »  CPLT 

IF    (NONSTK   .EQ.    0)    GO  TO  1 

KIM   =    LAND  -  NONSTK   +   1    +  LAST 
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IF    (KIM    . GT .    LAND)    KIM   =   KIM  -  LAND 
ICH   =   KIM  +   NONSTK   -  1 

C 

C  AGE   ZERO  OF   AREAS   SUBROUTINE   GIVES   CONDITIONS   AT   END   OF   YEAR  ZERO, 

C  SO   INCREASE  AGES  ONE  YEAR  TO  START  SIMULATION 

C 

1  DO  2    1=1, LAND 

AC  AGE ( I )    =   AC  AGE ( I )    +  1.0 

2  CONTINUE 

IF    (NONSTK    .EQ.    0)    GO   TO  4 

C 

C  SUPPRESS   AGE   INCREASE  FOR  NONSTOCKED  ACRES 

C 

DO  3  I=KIM, ICH 
AC  AGE (  I  )   =  0.0 

3  CONTINUE 

4  IF    (DEFOR   .EQ.   0.0)   GO  TO  15 

C 

C  DETERMINE   AREA  DEFORESTED  ANNUALLY 

C 

AKDX   =   LAND  -  NONSTK 

YRLOS  =   (AKDX  *  DEFOR)   +  YRLOS 

IF    (YRLOS    .LT.    1.0)    GO  TO  15 

C 

C  GENERATE   PSEUDORANDOM   NUMBER   FOR    AGE   OF    ACRE  DESTROYED 

C 

5  NDIV   =    (17.0   *   ANUL   +  3.01/128.0 
NULL  =  ANUL 

NULL   =    (17   »   NULL    +   3)    -   128   *  NDIV 
ANUL   =  NULL 

C 

C  CHECK  THAT  AGE   EXISTS  AND   IS  BETWEEN  ONE  AND  OLDEST  CURRENT  AGE 

C 

IF    (ANUL    .LE.    0.0)    GO   TO  5 

IF    (ANUL   .GT.    AC AGE ( KOUNT ) )   GO   TO  5 

DO   6   M=  1  ,  LAND 

K  AC  R   =  M 

IF    ( AC AGE ( M )    .EQ.   ANUL)    GO  TO  7 

6  CONTINUE 
GO  TO  5 

C 

C  SET  LOSS   TO  REDUCE  CURRENT   ALLOWABLE  CUT 

C 

7  LOSS  =   LOSS  +  1 
NONSTK   =   NONSTK   +  1 
YRLOS   =   YRLOS   -  1.0 

C 

C  SALVAGE  VOLUME    IF  NOT  LESS  THAN  BFSALV 

C  NO   SALVAGE   OR   CLEANUP    IF   AGE  LESS   THAN  AGEO. 

C 

IF    (NULL   .LE.   Nl)    GO  TO  9 

IF    (IYEAR   .EQ.    1)   MTHN  =  FMRCHD ( 1 ) 

NULL   =  NULL   +  1 

KULL   =   NULL   -  1 

IF    (KULL    .LT.    MTHN)    GO  TO  8 

SALVB  =  ANBDF ( NULL )    +   BDFC ( KULL ) 

IF    (SALVB   .LT.   BFSALV)    SCLOSS  =  SCLOSS  +  CLOSS 
IF    (SALVB   .LT.   BFSALV)    GO  TO  9 
VLBF   =   VLBF   +  SALVB 

RETTH   =   RETTH  +   SALVB   ♦    (PRIBDI IYRM)    »  BFPCT) 
GO   TO  9 

8  SALVB  =   ANBDF ( NULL ) 

IF    (SALVB   .LT.   BFSALV)    SCLOSS   =  SCLOSS   +  CLOSS 
IF    (SALVB   .LT.    BFSALV)   GO  TO  9 
VLBF   =   VLBF   +  SALVB 

RETTH   =   RETTH  +   SALVB   *    (PRIBDI IYRM)    *  BFPCT) 

C 

C  RENUMBER  ACRES   TO  PUT  ACRE   LOST   AT  END  OF  AGE  SEQUENCE 

C  WITH  AGE  ZERO 

C 

9  IF    (KACR    .NE.    KOUNT)    GO   TO  10 
LAST  =  LAST  +  1 

KOUNT  =  KOUNT  +  I 
AC  AGE ( LAST )  =  0.0 
GO  TO  15 

10  LUB  =   LAST   -  1 

IF    (KACR    .LT.    LAST)    GO   TO  13 
MNO  =   LAND  -  KACR 
DO    11  J=1,MN0 
JSUB  =  KACR  +  J 
1SUB  =    JSUB  -  1 
ACAGE(ISUB)    =  ACAGE(JSUB) 

11  CONTINUE 

AC AGE ( L  AND )   =  ACAGE(l) 

00   12  K=1,LUB 

KAN  =  K  +  1 

AC AGE ( K )   =  AC AG  E ( K AN ) 

12  CONTINUE 

ACAGE ( L AST )    =  0.0 
GO  TO  15 

13  00    14   M  =  K ACR , LUB 
MOL    =   M  +  1 

ACAGE ( M )   =   ACAGE ( MOL ) 

14  CONTINUE 

ACAGE ( LAST )   =  0.0 

IF    (YRLOS   .GE.    1.0)   GO  TO  5 

C 

C  PREPARE   SUBTOTALS   FOR  CURRENT  YEAR  AND  CHECK   THAT  NO  ACRE 

C  IS  OLDER  THAN   179  YEARS 

C 

15  DO    16  K=l,180 
NACRE ( K )    =  0 

16  CONTINUE 

DO    19   K= 1 , LAND 

IF    I  AC AGE ( K )    .LE.    179.0)    GO   TO  18 
WRITE  (6,17) 

17  FORMAT    I  1H1  ,/////, 47X.38HY0U  WENT  BEYOND  AGE  LIMIT  OF   179  YEARS 
IYEAR   =  NOYRS 


GO   TO  42 

18  LM  =   AC  AGE ( K )    +  1.0 
NACRE ( LM )   =  NACRE ( LM )   +  1 

19  CONTINUE 

C 

C  DETERMINE   ALLOWABLE  CUT  BASED  ON  BD.FT.    STUMPAGE  PRICE 

C 

DO  20  J=l,10 
NSUB  =  J 

IF    (PRIBD(IYRM)    .LE.    PRIDIV(J))    GO   TO  21 

20  CONTINUE 

21  IALCUT   =   MALCUT ( NSUB )    -  LOSS 
CUTAGE   =    FMRCHD ( NSUB ) 

C 

C  COMPUTE   THINNINGS  FOR  ANNUAL  CUT 

C 

MXY  =  0 

MAC   =  CUTAGE 

DO  24   I =N 1 ,MAC, JCYCL 

VBTH  =  0.0 

VCTH   =  0.0 

IF    (I    .GE.   MAC)   GO  TO  25 
MR   =    I    +  1 

IF    (BDFC(I)    .LT.    COMBF)    GO   TO  22 
VLBF 1   =   NACRE ( MR )    «  BDFC(I) 
VLBF   =   VLBF   +  VLBF1 
VBTH  =   VLBF 1 

RETTH   =   RETTH  +   VBTH   »    (PRIBDI IYRM)    *  BFPCT) 
MXY   =   MXY  +  1 
GO   TO  24 

22  IF    (CFMC(I)    .LT.   COMCU)    GO  TO  23 
VLCU1   =   NACRE  I  MR )    *  CFMC(I) 
VLCU  =  VLCU  +  VLCU1 

VCTH  =  VLCU1 

RETTH   =   RETTH   +   VCTH/100.0   *    (PRICF(IYRM)    *  CFPCT) 
MXY   =   MXY  +  1 
GO   TO  24 

23  MXY   =   MXY  +  1 

SCTHN   =   NACRE  I  MR )    *   CTHN   +  SCTHN 

24  CONTINUE 

25  MTHN   =   Nl   +   MXY   »  JCYCL 

IF    ( IALCUT    .LE.    0)    GO  TO  33 

C 

C  COMPUTE  VOLUME  OF  ACRES  HARVESTED 

C 

KYR    =    IYEAR   +   1    +   F I NL 

SHWD  =  SHELT  +   (SHELT  *  GROW) 

ISHWD  =   SHWD  +  0.5 

DO   32    1=1, IALCUT 

VBHV   =  0.0 

VCHV   =  0.0 

IF    (LAST   .LT.   LAND)   GO  TO  26 
LAST   =  0 

26  LAST   =   LAST   +  1 

IF    ( ACAGE ( L AST )    .GE.   CUTAGE)   GO  TO  27 
LAST   =   LAST  -  1 
GO  TO  33 

27  M  =  ACAGE ( L AST ) 
K  =  M  +  1 

KOUNT  =  KOUNT  +  1 
ISAFE   =   LANO  +  1 

IF    (KOUNT   .GE.    ISAFE)    KOUNT  =  1 

IF    I M   .LT.   MTHN)   GO  TO  28 

VLBF2  =   ANBDF ( K )   +  BDFC(M)   -  SHELT 

IF    ( VLBF2   .LT.   COMBF)   GO  TO  29 

VLBF   =   VLBF   +  VLBF2 

VBHV  =  VLBF2 

RETHV  =  RETHV  +  VBHV  *  PRIBDI IYRM I 
I VAR I  4 , KYR )  =  I VAR I  4, KYR )  +  ISHWD 
VAR ( 3 , KYR )   =  VAR ( 3 , KYR )   +    (SHWD  *  PR  I BD ( KYR ) ) 

C 

C  INCLUDE   FOLLOWING  STATEMENT   IF  REGENERATION  COSTS,    INCLUDING  FINAL 

C  CUT,   NOT  ALL   INCLUDED   IN  BFCST  AS  ASSESSMENT  AGAINST  VOLUME  OF 

C  CURRENT  HARVEST. 

C  VAR I  7 , KYR )   =  VARI7.KYR)    +   SHWD  *   (BFCST  +  BFCST  *  RATE  *  FIND 

C 

GO  TO  31 

28  VLBF2  =   ANBDF I K )  -  SHELT 

IF    (VLBF2   .LT.   COMBF)    GO  TO  29 
VLBF  =  VLBF  +  VLBF2 
VBHV  =  VLBF2 

RETHV  =  RETHV  +  VBHV  *  PRIBDI IYRM) 
IVARI4.KYR)  =  IVAR(4,KYR)  +  ISHWD 
VAR ( 3 , KYR )   =  VARI3.KYR)    +   (SHWD  *  PRIBDI KYR ) ) 

C 

C  INCLUDE   STATEMENT  FOR  VARI7.KYR)    HERE    IF    INCLUDED  JUST  PRIOR  TO 

C  STATEMENT  28,  ABOVE. 

C 

GO  TO  31 

29  IF    (M   .LT.   MTHN)   GO  TO  30 
VLCU2  =  ANCUV(K)   +  CFMC(M) 

IF   (VLCU2  .LT.  COMCU)   GO  TO  32 
VLCU  =  VLCU  +  VLCU2 
VCHV  =  VLCU2 

RETHV  =  RETHV  +  VCHV/100.0  *  PRICF(IYRM) 
GO  TO  31 

30  VLCU2  =  ANCUVIK) 

IF    (VLCU2   .LT.    COMCU)   GO  TO  32 
VLCU  =  VLCU  +  VLCU2 
VCHV  =  VLCU2 

RETHV  =  RETHV  +  VCHV/100.0  *  PRICFIIYRM) 

31  ACAGE I  LAST )   =  0.0 

32  CONTINUE 

C 

C  COMPUTE  GROWING  STOCK  VOLUME.   USE  CU.FT.    IF  VOLUME   IS  LESS 

C  THAN  BFMRCH  M  BD.  FT. 

C 

33  DO  37   1=1, LAND 
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IF    ( AC AGE ( I )    .IT.    AGMRCH)    GO  TO  37 
t  AG   =   AC AGE ( I )    +  1.0 
IBG   =    I  AG  -  1 

IF    IIRG    .LT.    MTHN)    GO   TO  34 
GBL1   =   ANBDF (  I  AG )    +  BDFC(IBG) 
IF    (GBL1    .LT.    BFMRCH)    GO   TO  35 
GVLBF   =   GVLBF   +  GBL1 
GO   TO  37 

34  GBL1   =   ANBDF ( I  AG ) 

IF    (GBL1    .LT.    BFMRCH)    GO   TO  35 
GVLBF   =   GVLBF   +  GBL1 
GO   TO  37 

35  IF    (IBG    .LT.    MTHN)    GO  TO  36 
GCL1   =   ANCUV(IAG)    +  CFMC(IBG) 
GVLCU   =   GVLCU   +  GCL1 

GO   TO  37 

36  GCL1    =   ANCUVI I  AG ) 
GVLCU   =   GVLCU   +  GCL1 

37  CONTINUE 

PREPARE   FOR   NEW   TOTALS   AND  SUBTOTALS 

DO   38  K=l,130 

38  NACRE ( K )    =  0 
DO    39  1=1,18 

39  ISUM(I)    =  0 

DO  40   K= 1 , L  AND 

LM   =   AC AGE ( K )    +  1.0 

NACRE(LM)    =   NACRE(LM)    +  1 

40  CONTINUE 

COMPUTE   TOTAL   ACREAGE    BY    10-YEAR   AGE  CLASSES 

DO   41  1=1,18 

DO   41  J=l,10 

NS   =    10   *    ( I   -    1 )    +  J 

ISUMl I )    =    ISUMI I )    +  NACRE(NS) 

41  CONTINUE 

COMPUTE  VOLUMES  AND  VALUES  AT  END . OF  CURRENT  YEAR  FOR  TRANSFER 
TO  ANUAL 


RETRN  =  RETTH  ^ 
CSTAC  =  LAND  * 
CSTVL  =  CUCST  ' 
TCOST  =  CSTAC  ^ 
ANNET  =  RETRN  - 
GSVALB   =  GVLBF 


RETHV 
ACCST   +   SCPLT  • 

(VLCU/100.0)-  • 

CSTVL 

TCOST 
*  (PRIBD(IYRM) 


SCTHN  +  SCLOSS 
BFCST   *  VLBF 


GSVALC   =    (GVLCU/100.0 )    *  (PRICF(IYRM) 


INCREASE   COSTS   ANNUALLY,    IF  DESIRED 


+    (ACCST   *  RATE) 
+    (BFCST   *  RATE) 
+    (CLOSS   *  RATE) 
(CPLT   *  RATE) 
(CTHN   *  RATE) 
+    (CUCST  *  RATE) 


ACCST  =  ACCST 
BFCST  =  BFCST 
CLOSS  =  CLOSS 
CPLT  =  CPLT  + 
CTHN  =  CTHN  + 
CUCST  =  CUCST 
RETURN 
ENO 

SUBROUTINE  0UTPUT1 

COMMON   ACCST, AGEO, AGMRCH, ANBDF ( 1 8 1 ) , ANC UV ( 1 8 1 ) , ANNE T , B A TC H ( 3 ) , BDFC 
1(180) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 180) , C ST AC , CS T V L , CUC ST , CUT AGE 
2,0BH0,DEN0,DESCR( 5 ),DLEV,FMRCHD(10) , GMN AM  I  3 ) , GS V AL B , GS V ALC , G VLBF , 
3GVLCU, I  ACRE! 1 80 ) , I ALC UT , I G AME , I SUM( 18) , I TE ST , I VAR ( 2 6,  1 50 ) , I  YEAR, 
4K0LI6) , LAND, LAST, MALCUT1 10) .MOLD, NACRE ( 1 80 ) , NG AME , NKOL S , N 1 , NONS TK , 
5N0YRS,PRET,PRIBD( 150), PRICFI 150), PR IDIV( 10), RETRN, ROTA, RATE, GROW, 
6SITE,SUMM(6,25, 10) , TCOST, SHELT , YRLOS , F I NL, CLOSS, CTHN, CPLT, CFPCT, 
7VAR (14,150) ,VLBF,VLCU, AC  AGE ( 1000) , C YCL , KOUNT , D EFOR , ANUL , M I  X , I PLNT , 
8R I  NT, THIN, BFMRCH, BF S AL V , COMCU , COMBF , BFPC T , GNTR 

WRITE   TABLE  HEADINGS 

WRITE  (6,1) 

L   FORMAT    ( 1H1,////,46X,28HDISTRIBUTI0N   OF   ACRES   BY  AGE) 

WRITE    (6,2)    ( BATCH (11,1=1,3) 
I   FORMAT    (1H    ,45X,7HBATCH  ,3A8) 

WRITE    (6,3)    I  TEST 
5   FORMAT    (1H   , 45X , 4HTE ST , I  4) 

WRITE    16,4)    (GMNAMI I ) , 1=1, 3) 
\   FORMAT    (1H    ,45X,6HGAME  ,3A8) 

WRITE    (6,5)    (DESCRI  I  )  ,  1=1,5) 
>   FORMAT    (1H  ,45X,5A8) 

WRITE    (6,6)  IYEAR 
i   FORMAT    ( 1H    ,45X,16HYEAR   WITHIN   GAME, 16,//) 

WRITE  (6,7) 
7   FORMAT    (1H    , 55X , 9HAGE ( YE AR ) I 

WRITE  (6,81 

3  FORMAT  (1H  ,  4X  ,  1 1 HAGE  (  DEC  ADE  )  ,  7X  ,  1H0  ,  7X  ,  1H1  ,  7X  ,  1H2 ,  7X  ,  1H3  ,  7X  ,  1H4  ,  7 
IX, 1H5.7X, 1H6,7X,1H7,7X, 1H8.7X, 1H9 , 1  OX , 5HT0TAL , / / ) 

-YEAR   AGE  CLASS  AND   TOTALS  OF 


DO    10  J=l,18 

K  =  J  -  1 

NN   =    10   *   K   +  1 

WRITE  (6,9)  K, NACRE (NN) , NACRE (NN+1 ) , NACRE ( NN+2 ) , NACRE ( NN+3 ) , NACRE ( 
1NN+4) , NACRE ( NN+5 ) , NACRE (NN+6) ,NACRE (NN+7) ,NACRE(NN+8) ,NACRE( NN  +  9) , 
2  I  SUM ( J ) 

9   FORMAT    (1H    ,  1 1 1 , 5 X , 10 1  8 , I  1 5 , / ) 

10  CONTINUE 
RETURN 
END 

SUBROUTINE  ANUAL 

COMMON  ACCST, AGEO, AGMRCH, ANBDF ( 1 8 1 ) , ANC UV ( 1 8 1 ) , ANNE T , 5 ATC H ( 3 ) , BDFC 
1  I  180) ,BDFO(  180) ,BFCST,CFMC(  180)  ,CFMO( 1 80 ) , C ST AC , C ST VL , CUC ST , CUT  AGE 
2, DBHO.DENO, DESCRI 5 ) , DL E V , FMRCHD ( 1 0 ) , GMN AM ( 3 ) , GSVALB .GSVALC , GVLBF , 


3GVLCU, I  ACRE ( 180) , IALCUT, I  GAME, I  SUM ( 18 ) ,ITEST,IVAR(26,150),IYEAR, 
4K0LI6) , LAND, LAST, MA LCUTl 10), MOLD, N ACRE! 180) , NG AME , NKOLS , N 1 , NONS TK , 
5N0YRS,PRET,PRIBD( 150 ) , PR  I CF ( 1 50 ) , PR  I D I V ( 10 ) , RETRN , ROT  A , RATE , GROW , 
6SITE,SUMM(6,25, 10) .TCOST, SHELT, YRLOS, F I NL, CLOSS, CTHN, CPLT, CFPCT, 
7V AR ( 14, 150), VLBF, VLCU, AC  AG E( 1000 ),C YCL, KOUNT, DEFOR, ANUL, MIX,  I  PLNT, 
8R I  NT, TH I N.BFMRCH.BFSALV, COMCU, COMBF, BFPC T, GNTR 

CONVERT   VOLUME   AND   AREA   VALUES   TO  SUBSCRIPTED   VALUES   FOR   USE  BY 
REMAINING  SUBROUTINES 


K  =  IYEAR 

J  =    IYEAR  + 

IVARd.J)  = 

IVARI2.J)  = 

IVARI3.J)  = 
I BF  T   =  VLBF 

IVARI4.J)  = 

IVARI5.J)  = 

IVARI6.J)  = 

IVARI7.J)  = 

IVARI8.J)  = 

IVAR(9,J)  = 

IVAR(IO.J)  : 


0.5 


IALCUT 
CUT AGE 
VLCU 
+  0.5 

IVARI4.J)  + 
IVAR(5,K)  + 
IVARI6.K)  + 
GVLCU  +  0.5 
GVLBF  +  0.5 
IVAR(5,J)  + 

:     IVARI6.J)  H 


IBFT 

IVAR13, J ) 
I VAR ( 4, J ) 


IVARI7.J 
I VAR ( 8 , 


IVARI 11, J)  =  NONS  TK 
DO    1  1=1,14 


IVAR(N.J)    =  ISUM(I) 
IVARI26.JI    =  ISUM115) 


ENTER   MONEY   VALUES    IN   ARRAYS   FOR   REMAINING  SUBROUTINES 


PRICFI J) 
PR  I BD ( J ) 
VAR(3,J)  + 
VAR(4,K)  + 
CSTAC 

VARI6.K)  + 
VAR(7,J)  + 
VARI8.K)  + 
VAR I  5 , J )  + 
=  VAR( 10, K) 
=  VAR(3,J)  ■ 
=   VAR( 12, K) 


RETRN 
VAR ( 3 , J ) 

VARI5.J) 
CSTVL 
VAR( 7, J) 
VAR ( 7 , J ) 
+   VAR ( 9, J ) 
■  VAR(9,J) 
+   VARI 11, J) 


GSVALC    +  GSVALB 
VARI12.J)    +  VARI13.J) 


VAR I  1 , J ) 
VARI2.JI 
VAR (3, J) 
VAR(4,J) 
VAR ( 5 , J ) 
VAR I  6 , J ) 
VAR(7,J) 
VAR(8,J) 
VAR ( 9 , J  ) 
VAR ( 1 0 , J ) 
VARI 11, J) 
VARI 12, Jl 
VARI13.JI 
VARI 14, J) 
RETURN 
END 

SUBROUTINE  0UTPUT2 

COMMON   ACCST.AGEO, AGMRCH, ANBDF (181), ANCUVI 181), ANNET, BATCHI3), BDFC 
1(180) ,BDFO( 180) ,BFCST,CFMC(180) ,CFMO( 18 0 ) , C ST AC , C ST VL , CUC ST , CUT AGE 
2,DBH0,DEN0,DESCR ( 5) , DL EV , FMRCHD ( 1 0 ), GMN AM ( 3 ), GSVALB , GSVA LC , GVLBF , 
3GVLCU, I  ACRE! 180) , I ALCUT , I  GAME , I  SUM ( 18 ) , I  TEST , I  V AR ( 26 , 1 50 ) , I  YEAR, 
4K0LI6) , LAND, L AST, MALCUTI 10) .MOLD, NACRE ( 1 80 ) , NGAME , NKOL S , N 1 , NONSTK , 
5N0YRS,PRET, PRIBDI 150), PRICFI 150), PRIDIV(IO), RETRN, ROTA, RATE, GROW, 
6SI TE.SUMMI 6,25, 10) , TCOST, SHELT, YRLOS, F I NL , C LOSS , CTHN , CPLT , CFPCT , 
7VARU4, 150) , VLBF, VLCU, ACAGE(  1000)  , C YCL , KOUNT , DE FOR , ANUL , M I  X , I PLNT , 
8R I N T, TH I N,BFMRCH,BFSALV, COMCU, COMBF, BFPCT, GNTR 

PRINT   FIRST  PAGE 

N  =   NOYRS  +  1 
H  =  40 

1  DO   15  J=1,N 
LINE  =   J  -  1 

2  IF    (M    .LT.   40)    GO   TO  11 
M   =  0 

WRITE    (6,3)    (BATCH! I ), 1=1, 3 ) 

3  FORMAT    ( 1H1,//,46X,7HBATCH  ,3A8) 
WRITE    (6,4)  ITEST 

4  FORMAT    (1H   ,45X , 4HT EST, 14) 
WRITE    (6,5)    (GMNAMI I ), 1=1,3) 

5  FORMAT    < 1H    ,45X,6HGAME  ,3A8) 
WRITE    (6,6)    (DESCRI I ), 1=1, 5) 

6  FORMAT    (1H  ,45X,5A8) 
WRITE  (6,7) 

7  FORMAT    (1H  ,/) 
WRITE  (6,8) 

8  FORMAT    (1H    , 1 2 X , 9H ALLOWABL E , 5X , 7HCUTT I NG , 8X , 1  OH AC TU AL  CUT,10X,9HCU 
1MUL   CUT,10X,9HGRSTK   VOL , 12X , 9HT0T AL  VOL) 

WRITE  (6,9) 

9  FORMAT    (1H   , 2X , 4HYE AR , 9X , 3HCUT , 10X , 3HAGE , 7X , 6HCU . FT . , 5X , 3HMBF , 6X , 6 
1HCU.FT. , 5X,3HMBF,6X,6HCU.FT. , 5X , 3HMBF , 6X , 6HCU . FT . ,5X,3HMBF ) 

WRITE  (6,10) 

10  FORMAT    (1H   , 15X , 3H ( 1 ) , 1  OX , 3H ( 2 ) , 8X , 3H ( 3 ) , 7X , 3H ( 4 ) , 7X , 3H ( 5 ) , 7X , 3H ( 6 
1),7X,3H(7),7X,3H(8),7X,3H(9),6X,4H(10),//) 

11  WRITE    (6,12)    LINE, ( I VARI I , J) , 1=1, 10) 

12  FORMAT    (1H  ,16,112,113,112,19,3(111,19)) 
IF    (J    .LE.    1)   GO  TO  13 

M  =   M   +  1 

IF    (LL    .LT.    10)    GO   TO  14 

13  WRITE  (6,7) 
LL   =  0 

14  LL   =   LL   +  1 

15  CONTINUE 

PRINT   SECOND  PAGE 


DO  23  J=l, 
LINE  =  J  - 
IF    IM  .LT. 


40)    GO   TO  19 


WRITE  (6,3)  (BATCH! 

WRITE  (6,4)  ITEST 

WRITE  (6,5)  (GMNAMI 

WRITE  (6,6)  (DESCRI 

WRITE  (6,7) 

WRITE  (6,16) 


21 


16  FORMAT    (1H    ,  1 1 X , 3HNON , 46 X , 1 1H AG E  CLASSES) 
WRITE  (6,17) 

17  FORMAT    (1H    , 2X , 4HYE AR , 5X , 3HSTK , 1 1 4H        0-9      10-19      20-29      30-39  40 
l-/t9     50-59     60-69     70-79     80-89     90-99     100-109     110-119  120-129 
2   130-139  140-179) 

WRITE  (6,18) 

18  FORMAT    (1H    , 1 OX , 1 1 7H ( 1 1 )         (12)      (13)         114)         (15)         (16)  (17) 

1  (18)        119)        (20)        (21)  (22)  (23)  (24)  (25) 

2  (26)  ,//) 

19  WRITE    (6,20)   LINE,(IVAR(  I, J), 1  =  11, 26) 

20  FORMAT    (1H  ,16,217,16,817,18,419) 
IF    (J    .LE.    1)    GO   TO  21 

M   =   M  +  1 

IF    (LL    .LT.    10)   GO  TO  22 

21  WRITE  (6,7) 
LL  =  0 

22  LL   =   LL   +  1 

23  CONTINUE 

PRINT  THIRD  PAGE 

M  =  40 

DO  31  J=1,N 

LINE  =  J  -  1 

IF    (M   .LT.   40)    GO  TO  27 

M  =  0 

WRITE    (6,3)  (BATCH(I),I=1,3) 

WRITE    (6,4)    I  TEST 

WRITE    (6,5)    (GMNAM(I), 1=1,3) 

WRITE    (6,6)    (OESCRI I ) , 1=1,5) 

WRITE  (6,7) 

WRITE  (6,24) 

24  FORMAT    I 1H    , 18X , 14HSTUMPAGE    PR  I CE , 9 X , 1 5HSTUMP AGE    I NCOME , 1 3X , 1 OHARE 
1A  COSTS, 15X, 12HV0LUME  COSTS) 

WRITE  (6,25) 

25  FORMAT    (1H   , 2X , 4H YE AR , 9X , 1  OH 100  CU. FT. , 5X , 3HMBF , 6X , 6HANNUAL , 5X , 9HC 
1UMUL AT ED, 6X.6H ANNUAL, 5X , 9HCUMUL AT ED , 6X , 6HANNUAL , 5 X , 9HCUMUL ATED ) 

WRITE  (6,26) 

26  FORMAT    (1H   , 1 8X , 4H ( 2 7 ) , 8X , 4H ( 28 ) , 6X , 4H ( 2 9 ) , 9X , 4H ( 30 ) , 9X , 4H ( 3 1 ) , 9X , 
14H( 32) ,9X,4H(33) ,9X,4H(34) ,//) 

27  WRITE    (6,28)    L I NE , ( V AR ( I , J ) , 1  =  1 , 8  ) 

28  FORMAT    ( 1H   , 1 6 , F  1 6. 2 , F 1 1 . 2 , F 12 . 0 , F 1 1 , 0, 2 ( F 1 5 . 0 , F 1 1 . 0 )  ) 
IF    (J    .LE.    1)   GO  TO  29 

M   =   M  +  1 

IF    (LL    .LT.    10)    GO   TO  30 

29  WRITE  (6,7) 
LL  =  0 

30  LL  =   LL  +  1 

31  CONTINUE 


PRINT   FOURTH  PAGE 


DO  39  J=1,N 
LINE  =  J  -  1 
IF    (M   .LT.   40)    GO  TO  35 

M  =  0 

WRITE    (6,3)    (BATCH! I ) ,1=1,3) 

WRITE    (6,4)    I  TEST 

WRITE    (6,5)    ( GMNAM ( I ) , I =1 , 3 ) 

WRITE    (6,6)    (DESCRt I ) ,1=1,5) 

WRITE  (6,7) 

WRITE  (6,32) 

32  FORMAT    (1H   , 1 8X , 10HT0TAL  COST , 1 7X , 1 OHNE T    INCOME, 13X, 13HCURRENT  VAL 
1UE,9X,5HT0TAL) 

WRITE  (6,33) 

33  FORMAT    (1H   , 2X , 4HYE AR , 8X , 6HANNUAL , 5X , 9HCUMUL A T EO , 7X , 6H ANNUAL , 5X , 9H 
1CUMULATED,7X,13HGR0WING  STOCK, 7X, 9HNET  WORTH) 

WRITE  (6,34) 

34  FORMAT    (1H    , 1 5X , 4H ( 35 ) , 9X , 4H ( 36 ) , 10X , 4H ( 37 ) , 9X , 4H ( 38 ) , 1 3X , 4H ( 39 ) , 1 
14X,4H(40),//) 

35  WRITE    (6,36)    L I NE , ( V AR ( I , J ) , I =9 , 14 ) 

36  FORMAT    (1H   , 1 6 , F 1 4. 0 , F 12 . 0 , F 1 5 . 0 , F 1 2 . 0 , 2F 1 8 . 0 ) 
IF    (J    .LE.    1)   GO  TO  37 

M  =   M   +  1 

IF    (LL    .LT.    10)   GO  TO  38 

37  WRITE  (6,7) 

LL  =  0 

38  LL   =  LL   +  1 

39  CONTINUE 
RETURN 
END 

SUBROUTINE  WORTH 

COMMON  ACCST, AGEO, AGMRCH, ANBDFI 1 8 1 ) , ANCU V ( 1 8 1 ) , ANNE T , BATCH ( 3 ), BDFC 
11  180) ,BDFO( 180) ,BFCST,CFMC (180) ,CFMO( 180 ) , C ST AC , C ST VL , CUC ST , CUT AGE 
2,DBHO,DENO,DESCR(5) , DL E V , F MRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V ALB , GS V ALC , GVLB F , 
3GVLCU, I  ACRE! 180) , IALCUT, I  GAME, I  SUM ( 1 8 ) , I TE S T , I  V AR ( 26 ,  1 50 ) , I  YE AR , 
4K0LI6), LAND, LAST, MALCUT (10), MOLD, N AC RE(180),NGAME,NKOLS,Nl,NONSTK, 
5N0YRS.PRET, PR  I  80 ( 1 50 ) , PR  I CF ( 1 50 ) , PR  I D I  V ( 10 ) , RE TRN , ROT A , RATE , GROW , 
6SITE,SUMM(6,2  5,10) , TCOS T , SHELT , YRLOS , F I NL , CLOSS , C THN , CPL T , CF PCT , 
7VAR(  14,150) ,VLBF, VLCU, AC AGE ( 1000) , C YC L , KOUNT , DEFOR , ANUL , M I  X , I PLNT , 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT,GNTR 

DIMENSION  CRATE! 20) , DISCI  20) ,D I  SGI  20) ,D IS  I (20) , PRE V( 20) ,PWTH(20) ,R 
1ATI0I20, 150) 

DO   1  1=1,20 

CRATE  I  I  )   =  0.0 

DISCI  I  )   =  0.0 

DISGI I  1    =  0.0 

DISI  I  I  )    =  0.0 

PREVI  I  )  =0.0 

1  PWTHI I )    =  0.0 
DO   2  1=1,20 
DO   2  J=l,150 

2  RATIO! I , J )    =  0.0 


DO   3  1=1,19 
K  —  I  +  1 

3  CRATE ( K )   =  CRATE ( I  )    +  0.005 

PREPARE  INTEREST  TABLE  FOR  PERIOD  NOYRS 
DO  4  J=l,20 

FACTR  =    1.0  +  CRATE ( J ) 
DO   4   K=l, NOYRS 
YRS   =  K 

PROD  =  ALOG ( FACTR )   *  YRS 
RATIO(J.K)    =  EXP(PROD) 

4  CONTINUE 

DISCOUNT  GROWING  STOCK  VALUE  AT  NOYRS 

DO  5  L=l,20 
KL   =  NOYRS  +  1 

DISG(L)    =  VAR(  13, KD/RATIOIL, NOYRS) 

5  CONTINUE 

DISCOUNT  ANNUAL  COSTS  AND  RETURNS 

DO  7  M=l,20 
PRESC  =  0.0 
PRESI  =  0.0 
SPRSC  =  0.0 
SPRSI  =  0.0 
DO  6  N= 1 , NOYRS 
I   =  N  +  1 

PRESC  =  VAR19.I )/RATIO(M,N) 
PRESI    =  VARI3.I )/RATIO(M,N) 
SPRSC  =  SPRSC  +  PRESC 
SPRSI    =   SPRSI   +  PRESI 

6  CONTINUE 
DISKM)   =  SPRSI 
DISC(M)   =  SPRSC 

7  CONTINUE 

OBTAIN   PRESENT  WORTH  AT  EACH  RATE 


DISGI IJ) 
DISC!  IJ) 
*  100.0 


V AR I  1 3 i 


PREPARE  ARRAY  OF  ALTERNATIVE  RATES 
CRATE  I  1 )    =  0.010 


DO  8  IJ=1,20 
PREV!  IJ)   =  DISI I IJ)  * 
PWTHIIJ)    =   PREV(IJ)  - 
CRATE ( I  J)   =  CRATE ( I  J ) 

8  CONTINUE 
WRITE  (6,9) 

9  FORMAT    (1H1,////,53X,29HPRESENT  WORTH  AND  RATE  EARNED) 
WRITE    (6,10)  (BATCH(I),I=1,3) 

10  FORMAT    (1H   ,52X,7HBATCH  ,3A8) 
WRITE    (6,11)    I  TEST 

11  FORMAT    (1H   ,52X,4HTEST, 14) 
WRITE    (6,12)    (GMNAM! I ) , 1=1,3) 

12  FORMAT    (1H   ,52X,6HGAME  ,3A8) 
WRITE    (6,13)    (DESCRI I  1,1  =  1,5) 

13  FORMAT    (1H  ,52X,5A8) 
WRITE    (6,14)  NOYRS 

14  FORMAT    (1H   ,  52X  ,  15HYE ARS   IN  PER  I OD , I  5 , / / ) 
WRITE    (6,15)  VAR(13,1) 

15  FORMAT    (1H   , 1 IX , 33HVALUE  OF   INITIAL  GROWING  STOCK — $, F10. 
WRITE  (6,16) 

16  FORMAT    (1H   , 57X ,  38HVALUES  DISCOUNTED  TO  PRESENT  (DOLLARS),/) 
WRITE  (6,17) 

17  FORMAT    (1H   , 1 IX , 8HC OMPOUND , 1 5X , 6HFUTURE , 34X , 5HS TOCK , 36X , 3HNET ) 
WRITE  (6,18) 

18  FORMAT    (1H   , 13X , 4HR ATE , 1 6X , 7HGR0W I NG , 1 5X , 3HALL , 1 7X , 4HPLUS , 16X ,  3HAL 
1L, 15X.7HPRESENT) 

WRITE  (6,19) 

19  FORMAT    (1H   , 1 1 X , 9H (PERCENT), 14X,5HST0CK,14X,7HINC0MES,13X,7HI NCOME 
1S,14X,5HC0STS,15X,5HW0RTH,/) 

DO  21  1=1,20 

WRITE  (6,20)CRATE(I),DISG(I),DISI(I),PREV(I),DISC(I),PWTH(I) 

20  FORMAT    (1H    ,12X,F5.1,12X,5(F10. 2,10X1,/) 

21  CONTINUE 
RETURN 
END 

SUBROUTINE  SUMRY 

COMMON  ACCST, AGEO, AGMRCH, ANBDF! 181) , ANCUVI 181 ) , ANNE T , BATCH ( 3)  ,BDFC 
1( 180) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 1 80  )  , C ST AC , C ST VL , CUCST , CUTAGE 
2,DBH0,DEN0,DESCR(5) , OLE  V , FMRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V AL B , GSVALC , GVLBF , 
3GVLCU, I  ACRE (180), IALCUT, I  GAME, I  SUM (  1 8 ) , I TE ST , I VAR ( 26 , 1 50 )  ,  I  YEAR , 
4K0L (6) .LAND, LAST, MALCUT ( 10) , MOLD, NACRE (  180) , NGAME , NKOLS , N 1 ,  NONSTK, 
5N0YRS,PRET,PRIBD( 150) , PR  I C F ( 1 50 ) , PR  I  0 1 V ( 1 0 ), RE TRN , ROTA , RATE  ,  GROW , 
6SITE,SUMM(6,25, 10 ), TCOST , SHELT , YRLOS , F I NL , CLOSS  ,  C THN , CPLT , CFPC T , 
7VAR ( 14,  150) , VLBF, VLCU, AC AGE ( 1000) , CYCL , KOUNT , DEFOR , ANUL  ,  M I X ,  I  PLNT , 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT,GNTR 

CONVERT   IVAR(I.J)    AND  VAR(I.J)    TO  SUMM(I,J,K)    AT   END  OF  EACH  GAME 

LIM  =   10  +  NOYRS/10 
DO  4  I=1,NK0LS 
DO  4  J=1,LIM 
K  =   KOL! I ) 

IF    (J    .GT.    10)    GO  TO  1 
JJ  =   J   +  1 
GO  TO  2 

1  JJ   =    10  *    (J  -    10)    +  1 

2  IF    (K   .GT.   26)   GO  TO  3 
SUMM( I, J, IGAME)    =  IVAR(K,JJ) 
GO  TO  4 

3  K  =   K  -  26 

SUMMI I, J, IGAME)   =  VAR(K,JJ) 

4  CONTINUE 

WRITE   SUMMARY  TABLES    IF   ALL  GAMES  FINISHED 
IF    (IGAME   .LT.   NGAME)   GO  TO  17 


-  22  - 


WRITE   PAGE   HEADINGS.    SEPARATE   PAGE   FOR   EACH  COLUMN  OF  OUTPUT2 
IDENTIFIED   IN  INPUT1 

DO    16  I=l,NKOLS 
WRITE  (6,5) 

5  FORMAT    (  1H1,////,46X,26HC0MPARIS0N  OF  ALTERNATIVES) 
WRITE    (6,6)    (BATCH! I ), 1=1,3) 

6  FORMAT    (1H   ,45X,7HBATCH  ,3A8) 
WRITE    (6,7)    I  TEST 

7  FORMAT    (1H   , 45X  , 4HT EST,  14) 
WRITE    (6,8)    (DESCRI I), 1=1, 5] 

8  FORMAT    (1H  ,45X,5A8) 
K  =   KOL ( I ) 

WRITE   (6,9)  K 

9  FORMAT    (  1-H   ,45X , 8HCOLUMN  ,13,///) 
WRITE  (6,10) 

10  FORMAT  (1H  , 5X , 4HYEAR , 6X , 6HGAME  1,6X,6HGAME  2 , 6X , 6HGAME  3,6X,6HGAM 
IE  4,6X,6HGAME  5,6X,6HGAME  6,6X,6HGAME  7,6X,6HGAME  8,6X,6HGAME  9,6X 
2 , 7HGAME  10,//) 


WRITE  SUMMI I,J,K)  FOR  EACH  OF  FIRST  10  YEARS  AND  FOR  END  OF 
EACH  DECADE 


DO   16  J=l,25 
IF    (J    .GT.  10) 
JJ  =  J 


GO  TO  11 


GO   TO  12 

11  JJ   =   10  *    (J  -  10) 

12  IF    IM    .LT.    5)    GO   TO  14 
WRITE  (6,13) 

13  FORMAT    (1H  ,//) 
M  =  0 

14  WRITE    (6,15)    J  J, I SUMMI I , J,L) ,L= 

15  FORMAT    (1H   , 1 9 , F 1 1 . 0 , 9F 1 2 . 0 ) 

16  M  =  M+l 

17  RETURN 
END 
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APPENDIX  2 
Output  of  Test  Problem 


YIELDS  PER  ACRE  OF  MANAGED,   EVEN-AGED  STANDS  OF   PONDEROSA  PINE   IN  THE  BLACK  HILLS 


SITE 

INDEX  60, 

20-YEAR 

CUTTING  CYCLE,  DENSITY 

LEVEL 

100 

STAND  BEFORE  AND 

AFTER  THINNING 

PERIODIC  CUT 

AND  MORTALITY 

STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
MBF 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
MBF 

30 
30 

1000 
496 

110 
73 

4.5 
5.2 

20 
20 

818 
534 

201 
201 

0. 
0. 

504 

37 

284 

0 

0 

40 

496 

104 

6.2 

28 

1146 

684 

0. 

50 
50 

496 
286 

133 
86 

7.0 
7.4 

35 
35 

1866 
1209 

1381 
954 

0. 
0. 

210 

47 

657 

427 

0 

60 

286 

108 

8.3 

41 

1806 

1556 

1.790 

70 
70 

286 
200 

129 
99 

9.1 
9.5 

47 
47 

2501 
1930 

2241 
1756 

3.880 
3.510 

86 

30 

571 

485 

.370 

80 

200 

116 

10.3 

52 

2515 

2341 

5.990 

90 
90 

200 
141 

132 
100 

11.0 
11.4 

57 
57 

3144 
2407 

2940 
2255 

8.990 
7.560 

59 

32 

737 

685 

1.430 

100 

141 

115 

12.2 

60 

2902 

2730 

10.590 

110 
110 

141 
102 

130 
100 

13.0 
13.4 

63 
63 

3457 
2701 

3263 
2555 

13.650 
11.020 

39 

30 

756 

708 

2.630 

120 

102 

114 

14.3 

66 

3216 

3055 

13.960 

130 
130 

102 
76 

127 
100 

15.1 
15.5 

69 
69 

3741 
2998 

3565 
2863 

16.980 
13.910 

26 

27 

743 

702 

3.070 

140 

76 

112 

16.4 

71 

3445 

3297 

16.740 

150 

76 

123 

17.2 

73 

3887 

3723 

19.510 
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GROWING  STOCK  OF   MANAGED  BLACK  HILLS  PONDEROSA  PINE 
SITE   INDEX  60,     20-YEAR  CUTTING  CYCLE 
DENSITY  LEVEL-  120  AND  100 

VOLUMES  PRESENT  PER  ACRE  AT  END  OF  EACH  YEAR 
MERCHANTABLE  CUBIC  FEET 


DECADE 

0 

1 

2 

3 

5 

6 

7 

8 

9 

0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3 

201.0 

249.3 

297.6 

345.9 

394.2 

442.5 

490.8 

539.1 

587.4 

635.7 

4 

684.0 

753.7 

823.4 

893.1 

962.8 

1032.5 

1102.2 

1171.9 

1241.6 

1311.3 

5 

954.0 

1014.2 

1074.4 

1134.6 

1194.8 

1255.0 

1315.2 

1375.4 

1435.6 

1495.8 

6 

1556.0 

1624.5 

1693.0 

1761.5 

1830.0 

1898.5 

1967.0 

2035.5 

2104.0 

2172.5 

7 

1756.0 

1814.5 

1873.0 

1931.5 

1990.0 

2048.5 

2107.0 

2165.5 

2224.0 

2282.5 

8 

2341.0 

2400.9 

2460.8 

2520.7 

2580.6 

2640.5 

2700.4 

2760.3 

2820.2 

2880.1 

9 

2255.0 

2302.5 

2350.0 

2397.5 

2445.0 

2492.5 

2540.0 

2587.5 

2635.0 

2682.5 

10 

2730.0 

2783.3 

2836.6 

2889.9 

2943.2 

2996.5 

3049.8 

3103.1 

3156.4 

3209.7 

11 

2555.0 

2605.0 

2655.0 

2705.0 

2755.0 

2805.0 

2855.0 

2905.0 

2955.0 

3005.0 

12 

3055.0 

3106.0 

3157.0 

3208.0 

3259.0 

3310.0 

3361.0 

3412.0 

3463.0 

3514.0 

13 

2863.0 

2906.4 

2949.8 

2993.2 

3036.6 

3080.0 

3123.4 

3166.8 

3210.2 

3253.6 

14 

3297.0 

3339.6 

3382.2 

3424.8 

3467.4 

3510.0 

3552.6 

3595.2 

3637.8 

3680.4 

15 

3723.0 

-  23  - 
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THOUSANDS   Of-   BOARD  FEET 


0. 
0. 

1.790 
3.510 
5.990 
7.560 
10.590 
11.020 
13.960 
13.910 
16.740 
19.510 


0. 

.  1  79 
1.999 
3.758 
6.290 
7.863 
10.896 
11.314 
14.262 
14.193 
17.017 


0. 

.358 
2.208 
4.006 
6.590 
8.  166 
11.202 
11.608 
14.564 
14.476 
17.294 


.537 
2.417 
4.254 
6.890 
8.469 
11.508 
11.902 
14.866 
14.759 
17.571 


0. 


.716 
2.626 
4.502 
7.190 
8.772 
11.814 
12.196 
15. 168 
15.042 
17.848 


0. 
0. 

0. 
0. 

.895 
2.835 
4.750 
7.490 
9.075 
12. 120 
12.490 
15.470 
15.325 
16.125 


7.790 
9.378 
12.426 
12.784 
15.772 
15.608 
18.402 


0. 
0. 
0. 
0. 
0. 

1.253 
3.253 
5.  246 
8.090 
9.681 
12.732 
13.078 
16.074 
15.891 
18.679 


0. 
0. 
0. 
0. 

0. 

1.432 
3.462 
5.  494 
8.  390 
9.984 
13.038 
13.372 
16.376 
16. 174 
18.956 


0. 
0. 
0. 

1.611 
3.671 
5.742 
8.690 
10.287 
13.344 
13.666 
16.578 
16.457 
19.233 
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ALTERNATIVES    FOR    THIS  GAME 
BATCH      TEST  PROBLEM 
TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED, THINNED   AGE  30 


CRITICAL  PRICES 
ALLOWABLE  CUT 
MINIMUM   CUTTING  AGE 


NUMBER   OF   YEARS   PER  GAME 


ACRES    IN  WORKING  CIRCLE  915 

MINIMUM  VALUES   FOR   INCLUSION   IN  TOTALS 

AGE.    FOR   GROWING   STOCK  40 

M  BD.    FT.,    FOR  GROWING  STOCK  1.5 

CU.    FT.,    FOR   COMMERCIAL   CUT  400 

M   BO.    FT.,    FOR   COMMERCIAL   CUT  3.0 

M   BD.    FT. ,    FOR   SALVAGE  1.5 


COSTS    IN   FIRST   YEAR   OF  GAME 

PER  ACRE  (ANNUAL) 

PER    100   CU.    FT.  HARVESTED 

PER   M   BD.  FT. 

THIN   ONE  ACRE 

PLANT   ONE  ACRE 

CLEANUP   OF    ONE  ACRE 
RATE   OF    INCREASE    IN  COSTS 


.05 
1.56 
25.00 
30.00 
25.00 

.01 


ACRES    PLANTED  ANNUALLY 
PERCENT   OF   ACRES   LOST  ANNUALLY 
M  BD.    FT.    IN  SHELTERWOOD 
REGENERATION  PERIOD 


.040 
4.0 
10.0 


RELATIVE   VALUE   OF    INTERMEDIATE  CUTS 

STUMPAGE  PRICE,  CU.  FT.  1.00 
STUMPAGE   PRICE,    BD.    FT.  .85 


PSEUDORANDOM   NUMBER  GENERATOR 


21.0 
2222.0 


PAGE  TYPE  4 


INITIAL    DISTRIBUTION   OF   ACRES   BY  AGE 
BATCH      TEST  PROBLEM 
TEST  1 

GAME      EOUAL    AREAS   CUT  ANNUALLY 
MANAGED, THINNED  AGE  30  ' 
YEAR   WITHIN   GAME  0 


AGE ( DEC ADE ) 
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PAGE  TYPE  4  (continued) 


DISTRIBUTION   OF   ACRES   BY  AGE 
BATCH      TEST  PROBLEM 
TEST  1 

GAME      EQUAL    AREAS   CUT  ANNUALLY 
MANAGED, THINNED   AGE  30 
YEAR   WITHIN   GAME  30 


AGE ( DECADE ) 


PAGE  TYPE   5  BATCH      TEST  PROBLEM 

TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED, THINNED   AGE  30 


CUTTING 
AGE 
(2) 


ACTUAL  CUT 
CU.FT.  MBF 
(3)  (4) 


CUMUL  CUT 
CU.FT.  MBF 
(5)  (6) 


TOTAL  VOL 
CU.FT.  MBF 
(9)  (101 


130 
130 
130 
130 
130 
130 
130 
130 
130 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
14719 
16135 


16135 
32270 
48405 
64540 
80675 
96810 
1  12945 
129080 
143799 
159934 


91 
162 
272 
363 
454 
545 
636 
727 
805 
897 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 


4439 
4439 
4444 
4444 
4443 
4447 
4447 
4446 
4468 
4467 


161243 
177378 
193513 
209648 
225783 
241918 
258053 
274188 
288907 
305042 


4530 
4621 
4716 
4807 
4897 
4992 
5083 
5173 
5273 
5364 


130 
130 
130 
130 
130 
130 
130 
130 
130 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
15450 


114 

127 
122 
114 
127 
122 
122 
128 
123 
119 


176069 
192204 
208339 
224474 
240609 
256744 
272879 
289014 
305149 
320599 


1011 
1138 
1260 
1374 
1501 
1623 
1745 
1873 
1996 
2115 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
144424 
144354 


4483 
4483 
4482 
4499 
4498 
4497 
4506 
4505 
4504 
4517 


321177 
337312 
353447 
369582 
385717 
401852 
417987 
434122 
449573 
464953 


5494 
5621 
5742 
5873 
5999 
6120 
6251 
6378 
6500 
6632 


130 
130 
130 
130 
130 
130 
130 
130 
130 
130 


16135 
16135 
15427 
16135 
16135 
16135 
16135 
16135 
15708 
16135 


123 
128 
128 
123 
117 
128 
123 
118 
129 
124 


336734 
352869 
368296 
384431 
400566 
416701 
432836 
448971 
464679 
480814 


2238 
2366 
2494 
2617 
2734 
2862 
2985 
3103 
3232 
3356 


144284 
144215 
144145 
144075 
144006 
143936 
143866 
143797 
144154 
144094 


4516 
4519 
4520 
4518 
4531 
4530 
4528 
4540 
4539 
4537 


481018 
497084 
512441 
528506 
544572 
560637 
576702 
592768 
608833 
624908 


6754 
6885 
7014 
7135 
7265 
7392 
7513 
7643 
7771 
7893 
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PAGE  TYPE  5  (continued)  batch    test  problem 

TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGEO, THINNED  AGE  30 


NON  AGE  CLASSES 


EAR 

STK 

0-9 

70-79 

100—109 

110-119 

120-129 

130-1 39 

140-17 

(11) 

<  12 ) 

2?i4? 

3?i  5? 

*?16? 

5(*i7? 

6(*i8? 

(  19 ) 

8  (*20 1 

(  21 ) 

(  22 ) 

(  23 ) 

(  24) 

( 25 ) 

( 26) 

70 

o 

, 

75 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

0 

Q 

2 

3 

75 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

o 

o 

3 

3 

75 

70 

70 

70 

70 

70 

70 

70 

70 

70 

69 

70 

70 

1 

o 

^| 

5 

70 

70 

70 

70 

70 

70 

70 

70 

70 

69 

70 

70 

75 

70 

70 

70 

70 

70 

70 

70 

70 

70 

68 

70 

70 

2 

- 

7 

n 

75 

70 

70 

70 

70 

70 

70 

70 

70 

70 

68 

70 

70 

2 

9 

75 

70 

70 

70 

70 

70 

70 

70 

70 

68 

70 

70 

2 

o 

9 

1 

75 

70 

70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

3 

0 

10 

74 

71 

70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

3 

o 

1 1 

72 

72 

70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

_ 

'J 

z? 

70 

69 

70 

70 

70 

70 

70 

70 

70 

fc« 

70 

7(1 

70 

70 

70 

70 

70 

70 

68 

70 

g 

n 

69 

75 

69 

69 

70 

70 

70 

70 

70 

70 

70 

AS 

70 

jj 

1  6 

69 

75 

69 

69 

70 

70 

70 

70 

70 

68 

70 

I  7 

l 

69 

75 

69 

69 

70 

70 

70 

70 

69 

70 

70 

68 

70 

jj 

69 

69 

69 

70 

69 

70 

68 

19 

69 

75 

69 

70 

69 

70 

70 

70 

69 

70 

70 

70 

68 

6 

o 

20 

70 

73 

70 

70 

69 

70 

70 

69 

70 

69 

70 

70 

68 

7 

0 

21 

70 

72 

71 

70 

69 

70 

70 

69 

70 

69 

70 

70 

68 

7 

0 

22 

70 

"71 

72 

70 

69 

70 

70 

69 

70 

69 

69 

70 

68 

8 

0 

23 

70 

71 

72 

70 

69 

70 

70 

69 

70 

69 

70 

69 

68 

a 

0 

24 

70 

70 

74 

69 

69 

70 

70 

69 

70 

69 

70 

69 

68 

8 

0 

25 

70 

69 

75 

69 

69 

70 

69 

69 

70 

69 

70 

69 

68 

9 

0 

26 

70 

69 

75 

69 

69 

70 

69 

70 

69 

69 

70 

69 

68 

9 

0 

27 

70 

69 

75 

69 

69 

70 

69 

70 

69 

69 

70 

69 

68 

9 

0 

28 

70 

69 

75 

69 

69 

70 

68 

70 

69 

69 

70 

69 

68 

10 

0 

29 

0 

70 

69 

75 

69 

70 

69 

68 

70 

69 

69 

70 

69 

70 

8 

0 

30 

0 

70 

70 

73 

70 

70 

69 

68 

70 

69 

70 

69 

69 

70 

8 

0 

BATCH     TEST  PROBLEM 
TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED, THINNED  AGE  30 


STUMPAGE  PRICE 

YEAR                    100  CU.FT.  MBF 

(27)  (28) 

0  2.50  14.50 

1  2.50  15.20 

2  2.50  17.80 

3  2.50  16.80 

4  2.50  13.40 

5  2.50  14.10 

6  2.50  17.40 

7  2.50  11.80 

8  2.50  11.10 

9  2.50  12.20 

10  2.50  12.90 

11  2.50  10.10 

12  2.50  8.30 

13  2.50  9.00 

14  2.50  10.90 

15  2.50  13.90 

16  2.50  13.10 

17  2.50  11.90 

18  2.50  12.70 

19  2.50  15.70 

20  2.50  13.60 

21  2.50  12.10 

22  2.50  15.20 

23  2.50  16.10 

24  2.50  16.70 

25  2.50  19.60 

26  2.50  18.50 

27  2.50  14.70 

28  2.50  15.50 

29  2.50  17.10 

30  2.50  13.00 


STUMPAGE   INCOME  AREA 


ANNUAL  CUMULATED  ANNUAL 

(29)  (30)  (31) 

0  0  0 

1784  1784  388 

2021  3805  392 

1880  5686  396 

1625  7310  400 

1688  8999  404 

1952  10950  408 

1482  12433  412 

1418  13851  384 

1325  15176  361 

1586  16762  424 

1566  18329  395 

1469  19798  433 

1512  21310  403 

1661  22971  407 

2192  25163  446 

2021  27184  416 

1862  29046  420 

2041  31087  460 

2346  33433  428 

2016  35449  433 

1904  37353  473 

2339  39692  441 

2471  42163  483 

2480  44642  419 

2726  47368  455 

2805  50173  498 

2235  52408  464 

2253  54661  468 

2617  57278  513 

2027  59305  511 


VOLUME  COSTS 


CUMULATED             ANNUAL  CUMULATED 

(32)  (33)  (34) 

0  0  0 

388  150  150 

780  151  301 

1176  151  452 

1575  155  607 

1979  156  763 

2387  158  921 

2799  160  1081 

3183  162  1242 

3543  140  1383 

3968  165  1548 

4363  145  1693 

4796  169  1862 

5199  171  2033 

5607  149  2182 

6053  175  2357 

6468  175  2533 

6888  169  2703 

7348  181  2883 

7776  182  3065 

8208  168  3233 

8682  187  3420 

9123  189  3608 

9606  190  3799 

10025  193  3992 

10479  174  4165 

10977  197  4362 

11440  199  4562 

11909  181  4742 

12421  203  4946 

12933  206  5152 
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PAGE  TYPE  5  (continued) 


BATCH     TEST  PROBLEM 
TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED.THINNED  AGE  30 


TOTAL  COST 


ANNUAL 
(35) 


CUMULATED 
(36) 


NET  INCOME 


CUMULATED 
(38) 


CURRENT  VALUE 
GROWING  STOCK 
(39) 


TOTAL 
NET  WORTH 
(40) 


538 
543 
547 
555 
560 
565 
572 
545 
501 
590 


538 
1081 
1628 
2182 
2742 
3308 
3880 
4425 
4926 
5516 


1247 
1477 
1334 
1070 
1128 
1386 
910 
873 
824 
997 


1247 
2724 
4058 
5128 
6256 
7643 
8553 
9426 
10250 
11247 


64099 
75570 
71139 
55955 
58989 
73641 
48657 
45464 
50507 
53551 


65346 
78294 
75197 
61083 
65246 
81284 
57210 
54890 
60757 
64798 


540 
602 
574 
557 
621 
592 
589 
640 
611 
600 


6056 
6658 
7232 
7789 
8410 
9002 
9591 
10231 
10841 
11442 


1026 
867 
938 
1104 
1571 
1429 
1273 
1401 
1736 
1415 


12273 
13140 
14078 
15182 
16753 
18182 
19455 
20856 
22592 
24007 


41104 
32952 
36007 
44593 
58001 
54315 
48919 
52431 
65832 
56440 


53377 
46092 
50085 
59775 
74755 
72497 
68374 
73287 
88424 
80447 


660 
630 
673 
611 
628 
695 
663 
649 
716 
717 


12102 
12731 
13405 
14016 
14644 
15339 
16002 
16651 
17367 
18084 


1244 
1709 
1798 
1868 
2098 
2110 
1572 
1604 
1901 
1310 


25252 
26960 
28758 
30626 
32724 
34834 
36406 
38010 
39911 
41221 


49565 
63512 
67506 
70106 
83350 
78246 
60917 
64611 
71762 
53038 


74816 
90472 
96264 
100732 
116074 
113080 
97324 
102621 
111673 
94259 


PAGE  TYPE  6 


PRESENT  WORTH  AND  RATE  EARNED 
BATCH     TEST  PROBLEM 
TEST  1 

GAME     EQUAL  AREAS  CUT  ANNUALLY 
MANAGED.THINNED  AGE  30 
YEARS   IN  PERIOD  30 


VALUE  OF   INITIAL  GROWING  STOCK- 


VALUES  DISCOUNTED  TO  PRESENT  (DOLLARS) 


COMPOUND 

RATE 
(PERCENT) 

1.0 

1-5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 

9.0 

9.5 
10.0 
10.5 


FUTURE 
GROWING 
STOCK 

39350.04 

33931.66 

29280.69 

25285.44 

21850.92 

18896.26 

16352.58 

14161.12 

12271.78 

10641.74 

9234.44 

8018.59 

6967.44 

6058.04 

5270.77 

4588.75 

3997.53 

3484.68 

3039.53 

2652.88 


50396.60 
46618.55 
43222.05 
40163.29 
37403.89 
34910.23 
32652.78 
30605.60 
28745.87 
27053.48 
25510.68 
24101.82 
22813.03 
21632.05 
20548.01 
19551.25 
18633. 19 
17786.21 
17003.50 
16278.99 


STOCK 
PLUS 
INCOMES 

89746.64 

.80550.21 

72502.74 

65448.73 

59254.81 

53806.49 

49005.36 

44766.72 

41017.65 

37695.22 

34745.12 

32120.41 

29780.47 

27690.10 

25818.77 

24140.00 

22630.72 

21270.89 

20043.02 

18931.87 


ALL 

COSTS 

15451.27 
14330.44 
13320.23 
12408.09 
11583.01 
10835.34 
10156.61 
9539.35 
8977.00 
8463.76 
7994.52 
7564.75 
7170.46 
6808.09 
6474.49 
6166.86 
5882.70 
5619.79 
5376.15 
5150.01 


NET 
PRESENT 
WORTH 

13303.36 

5227.76 

-1809.50 

-7951.37 

-13320.21 

-18020.87 

-22143.27 

-25764.64 

-28951.36 

-31760.55 

-34241.41 

-36436.36 

-38382.01 

-40110.01 

-41647.73 

-43018.87 

-44243.99 

-45340.91 

-46325. 14 

-47210.15 
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PAGE  TYPE  3  alternatives  for  this  game 

BATCH      TEST  PROBLEM 
TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED   AGE  30 


NUMBER   OF   YEARS   PER   GAME  30 


CRITICAL   PRICES  12.00  15.00         99.00  0.  0.  0.  0.  0.  0. 

ALLOWABLE  CUT  57  10  000000 

MINIMUM  CUTTING  AGE  130  130  120  0  0  0  0  0  0 


ACRES   IN  WORKING  CIRCLE  915 

MINIMUM  VALUES   FOR   INCLUSION   IN  TOTALS 

AGE,   FOR  GROWING   STOCK  AO 

M   BD.   FT.,   FOR  GROWING  STOCK  1.5 

CU.   FT.,    FOR  COMMERCIAL  CUT  400 

M   BD.    FT.,    FOR   COMMERCIAL    CUT  3.0 

M   BD.    FT. ,    FOR    SALVAGE  1.5 


COSTS   IN  FIRST  YEAR  OF  GAME 

PER   ACRE    (ANNUAL)  .20 

PER   100  CU.   FT.   HARVESTED  .05 

PER  M  BD.   FT.  1.56 

THIN  ONE   ACRE  25.00 

PLANT   ONE   ACRE  30.00 

CLEANUP  OF   ONE  ACRE  25.00 

RATE  OF   INCREASE   IN  COSTS  .01 


ACRES   PLANTED   ANNUALLY  1 
PERCENT  OF   ACRES   LOST   ANNUALLY  .040 

M   BD.    FT.    IN    SHE  L  TER WOOD  4.0 

REGENERATION   PERIOD  10.0 

PSEUDORANDOM  NUMBER  GENERATOR  21.0 


2222.0 


RELATIVE   VALUE   OF    INTERMEDIATE  CUTS 

STUMPAGE  PRICE,  CU.  FT.  1.00 
STUMPAGE   PRICE,    BD.   FT.  .85 


PAGE  TYPE  4  INITIAL   DISTRIBUTION  OF  ACRES  BY  AGE 

BATCH     TEST  PROBLEM 
TEST  1 

GAME     VARY  CUT  WITH  PRICE 
MANAGED, THINNED  AGE  30 
YEAR   WITHIN   GAME  0 


AGE( YEAR) 


DECADE) 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

TOTAL 

0 

12 

7 

7 

7 

7 

7 

7 

7 

75 

1 

7 

7 

7 

7 

7 

7 

7 

70 

2 

7 

7 

7 

7 

7 

7 

7 

70 

3 

7 

7 

7 

7 

7 

7 

7 

70 

4 

7 

7 

7 

7 

7 

7 

7 

70 

5 

7 

7 

7 

7 

7 

7 

7 

70 

6 

7 

7 

7 

7 

7 

7 

7 

70 

7 

7 

7 

7 

7 

7 

7 

7 

70 

8 

7 

7 

7 

7 

7 

7 

7 

70 

9 

7 

7 

7 

7 

7 

7 

7 

70 

10 

7 

7 

7 

7 

7 

7 

7 

70 

11 

7 

7 

7 

7 

7 

7 

7 

70 

12 

7 

7 

7 

7 

7 

7 

7 

70 

13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

17 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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PAGE  TYPE  4  (continued) 


DISTRIBUTION   OF    4CRES   BY  AGE 
B4TO      TEST  PROBLEM 
TEST  1 

GA"E     VARY  CUT  with  price 
MANAGED, THINNED   AGE  30 
YEAR   WITHIN  GAME  30 


AGE ( DECADE ) 


PAGE  TYPE  5  batch    test  problem 

TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED   AGE  30 


ALLOWABLE 

CUT 
(1) 


CUMUL  CUT 
CU.FT.  MBF 
(5)  16) 


GRSTK  VOL 
CU.FT.  MBF 
(7)  (8) 


TOTAL  VOL 
CU.FT.  MBF 
19)  (10) 


120 
120 
120 
130 
130 
120 
130 
130 
130 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
14719 
16135 


129 
128 
126 
0 
78 
129 
65 
65 
78 


16135 
32270 
48405 
64540 
80675 
96810 
112945 
129080 
143799 
159934 


129 
257 
383 
383 
461 
590 
655 
720 
798 
890 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 


4389 
4339 
4294 
4410 
4426 
4380 
4413 
4446 
4468 
4467 


161243 
177378 
193513 
209648 
225783 
241918 
258053 
274188 
288907 
305042 


4518 
4596 
4677 
4793 
4887 
4970 
5068 
5166 
5266 
5357 


130 
130 
130 
130 
130 
130 
130 
130 
120 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
15450 


103 
116 
112 

54 
125 
140 

87 
120 
166 
121 


176069 
192204 
208339 
224474 
240609 
256744 
272879 
289014 
305149 
320599 


993 
1109 
1221 
1275 
1400 
1540 
1627 
1747 
1913 
2034 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
144424 
144354 


4517 
4551 
4585 
4636 
4636 
4635 
4679 
4678 
4625 
4638 


321177 
337312 
353447 
369582 
385717 
401852 
417987 
434122 
449573 
464953 


5510 
5660 
5806 
5911 
6036 
6175 
6306 
6425 
6538 
6672 


130 
120 
120 
120 
120 
120 
130 
120 
120 
130 


16135 
16135 
15427 
16135 
16135 
16135 
16135 
16135 
15708 
16135 


115 
160 
159 
153 
156 
166 
111 
155 
178 


336734 
352869 
368296 
384431 
400566 
416701 
432836 
448971 
464679 


2149 
2309 
2468 
2621 
2777 
2943 
3054 
3209 
3387 
3430 


144284 
144215 
144145 
144075 
144006 
143936 
143866 
143797 
144154 
144094 


4637 
4588 
4537 
4484 
4446 
4393 
4392 
4353 
4302 
4401 


481018 
497084 
512441 
528506 
544572 
560637 
576702 
592768 
608833 
624908 


6786 
6897 
7005 
7105 
7223 
7336 
7446 
7562 
7689 
7831 
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PAGE  TYPE  5  (continued)  BATCH    TEST  problem 

TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED t THINNED  AGE  30 


NON  AGE  CLASSES 


BAR 

STK 

0-9 

10-19  20-29 

30-39 

40-49 

50-59     60-69     70-79     80-89  90-99 

100-109 

110-119     120-129  130-139 

140-1" 

(11) 

(12) 

(13)  (14) 

(15) 

(16) 

( 17) 

(  18) 

( 19) 

20 

(21) 

(22) 

(  23) 

( 24) 

( 25) 

(26 

0 

5 

75 

70  70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

0 

0 

1 

4 

78 

70  70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

67 

o 

0 

2 

3 

81 

70  70 

70 

70 

70 

70 

70 

70 

70 

70 

70 

64 

0 

0 

3 

84 

70  70 

70 

70 

70 

70 

70 

70 

70 

69 

70 

62 

0 

0 

4 

77 

70  70 

70 

70 

70 

70 

70 

70 

70 

69 

70 

69 

o 

0 

5 

I 

76 

70  70 

70 

70 

70 

70 

70 

70 

70 

69 

70 

70 

0 

0 

6 

I 

79 

70  70 

70 

70 

70 

70 

70 

70 

70 

68 

70 

68 

o 

0 

7 

77 

70  70 

70 

70 

70 

70 

70 

70 

70 

68 

70 

70 

o 

0 

8 

75 

70  70 

70 

70 

70 

70 

70 

70 

70 

68 

70 

70 

2 

0 

9 

75 

70  70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

3 

0 

10 

74 

71  70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

3 

0 

11 

I 

67 

75  70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

6 

0 

12 

61 

79  70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

8 

0 

13 

56 

82  70 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

10 

0 

1* 

I 

60 

76  69 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

13 

0 

15 

60 

76  69 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

13 

0 

16 

57 

79  69 

69 

70 

70 

70 

70 

70 

70 

70 

68 

70 

13 

0 

17 

1 

57 

77  69 

69 

70 

70 

70 

70 

69 

70 

70 

68 

70 

16 

0 

18 

59 

75  69 

69 

70 

70 

70 

70 

69 

70 

70 

68 

70 

16 

0 

19 

62 

75  69 

70 

69 

70 

70 

70 

69 

70 

70 

70 

68 

13 

0 

20 

63 

73  70 

70 

69 

70 

70 

69 

70 

69 

70 

70 

68 

14 

0 

21 

65 

67  74 

70 

69 

70 

70 

69 

70 

69 

70 

70 

68 

14 

0 

22 

70 

61  78 

70 

69 

70 

70 

69 

70 

69 

69 

70 

68 

12 

0 

23 

75 

56  81 

70 

69 

70 

70 

69 

70 

69 

70 

69 

68 

9 

0 

24 

80 

60  76 

69 

69 

70 

70 

69 

70 

69 

70 

69 

68 

6 

0 

25 

; 

83 

60  76 

69 

69 

70 

69 

69 

70 

69 

70 

69 

68 

4 

0 

26 

86 

57  79 

69 

69 

70 

69 

70 

69 

69 

70 

69 

68 

1 

0 

27 

88 

57  77 

69 

69 

70 

69 

70 

69 

69 

70 

69 

68 

1 

0 

28 

y 

91 

59  75 

69 

69 

70 

68 

70 

69 

69 

70 

69 

67 

0 

0 

29 

0 

91 

62  75 

69 

70 

69 

68 

70 

69 

69 

70 

69 

64 

0 

0 

30 

0 

85 

63  73 

70 

70 

69 

68 

70 

69 

70 

69 

69 

70 

0 

0 

BATCH  TEST 

PROBLEM 

TEST 

1 

GAME 

VARY  CUT  WITH 

PRICE 

MANAGED, THINNED  AGE 

30 

STUMPAGE 

PRICE 

STUMPAGE 

INCOME 

AREA  COSTS 

VOLUME 

COSTS 

YEAR 

100  CU.FT. 

MBF 

ANNUAL 

CUMULATED 

ANNUAL 

CUMULATED 

ANNUAL 

CUMULATED 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

( 34) 

0 

2.50 

14.50 

0 

0 

0 

0 

0 

0 

1 

2.50 

15.20 

2363 

2363 

388 

388 

209 

209 

2 

2.50 

17.80 

2682 

5044 

392 

780 

210 

419 

3 

2.50 

16.80 

2489 

7533 

396 

1176 

209 

628 

4 

2.50 

13.40 

403 

7936 

400 

1575 

8 

636 

5 

2.50 

14.10 

1501 

9438 

404 

1979 

135 

771 

6 

2.50 

17.40 

2614 

12052 

408 

2387 

220 

991 

7 

2.50 

11.80 

1169 

13221 

412 

2799 

116 

1107 

8 

2.50 

11.10 

1124 

14345 

384 

3183 

117 

1224 

9 

2.50 

12.20 

1325 

15670 

361 

3543 

140 

1364 

10 

2.50 

12.90 

1586 

17256 

424 

3968 

165 

1530 

11 

2.50 

10.10 

1462 

18718 

395 

4363 

100 

1629 

12 

2.50 

8.30 

1385 

20103 

433 

4796 

124 

1754 

13 

2.50 

9.00 

1424 

21527 

403 

5199 

126 

1880 

14 

2.50 

10.90 

991 

22518 

407 

5607 

105 

1985 

15 

2.50 

13.90 

2151 

24669 

446 

6053 

179 

2164 

16 

2.50 

13.10 

2255 

26924 

416 

6468 

181 

2344 

17 

2.50 

11.90 

1439 

28363 

420 

6888 

123 

2468 

18 

2.50 

12.70 

1945 

30308 

460 

7348 

186 

2653 

19 

2.50 

15.70 

3024 

33333 

428 

7776 

263 

2916 

20 

2.50 

13.60 

2039 

35372 

433 

8208 

171 

3087 

21 

2.50 

12.  10 

1802 

37174 

473 

8682 

190 

3278 

22 

2.50 

15.20 

2815 

39989 

441 

9123 

269 

3547 

23 

2.50 

16.10 

2962 

42952 

483 

9606 

270 

3816 

24 

2.50 

16.70 

2970 

45922 

419 

10025 

271 

4087 

25 

2.50 

19.60 

3484 

49406 

455 

10479 

250 

4338 

26 

2.50 

18.50 

3499 

52905 

498 

10977 

272 

4610 

27 

2.50 

14.70 

2049 

54954 

464 

11440 

195 

4804 

28 

2.50 

15.50 

2821 

57774 

468 

11909 

256 

5060 

29 

2.50 

17.10 

3470 

61245 

513 

12421 

274 

5334 

30 

2.50 

13.00 

978 

62223 

511 

12933 

38 

5372 
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PAGE  TYPE  5  (continued)  batch    test  problem 

TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED  AGE  30 


TOTAL  COST  NET   INCOME  CURRENT  VALUE  TOTAL 

YEAR  ANNUAL  CUMULATED  ANNUAL  CUMULATED  GROWING  STOCK  NET  WORTH 


PAGE  TYPE  6 


(35) 

(36) 

(371 

(38) 

(39) 

(40) 

0 

0 

0 

0 

0 

60992 

60992 

1 

2 
3 
4 
5 
6 
7 
8 
9 
10 

597 
602 
604 
408 
539 
628 
528 
501 
501 
590 

597 
1199 
1803 
2211 
2750 
3378 
3905 
4406 
4908 
5497 

1765 
2080 
1885 
-5 
963 
1986 
641 
623 
824 
997 

1765 
3845 
5730 
5725 
6688 
8674 
9316 
9938 
10762 
11759 

63417 
73946 
68856 
55558 
58777 
72580 
48313 
45464 
50507 
53551 

65182 
77792 
74586 
61283 
65465 
81255 
57628 
55403 
61269 
65310 

11 
12 
13 
14 
15 
16 
17 
18 
19 
20 

495 
557 
530 
512 
625 
596 
543 
645 
691 
603 

5993 
6550 
7079 
7592 
8216 
8813 
9356 
10001 
10692 
11296 

966 
828 
894 
479 
1526 
1659 
896 
1300 
2333 
1435 

12725 
13553 
14448 
14926 
16453 
18111 
19008 
20307 
22641 
24076 

41389 
33399 
36751 
45848 
59666 
55867 
50660 
54307 
67500 
57855 

54114 
46953 
51199 
60774 
76119 
73979 
69668 
74615 
90141 
81931 

21 
22 
23 
24 
25 
26 
27 
28 
29 
30 

664 
710 
753 
689 
705 
770 
658 
724 
787 
549 

11960 
12670 
13423 
14112 
14817 
15587 
16245 
16969 
17756 
18305 

1139 
2105 
2210 
2281 
2779 
2729 
1391 
2097 
2684 
429 

25214 
27320 
29529 
31810 
34589 
37318 
38709 
40806 
43489 
43918 

50797 
64429 
67750 
69601 
8  1844 
75995 
59188 
62090 
68202 
51552 

76011 
91749 
97279 
101411 
116433 
113313 
97897 
102896 
111691 
95470 

PRESENT  WORTH  AND  RATE  EARNED 
BATCH     TEST  PROBLEM 
TEST  1 

GAME     VARY  CUT  WITH  PRICE 
MANAGED, THINNED  AGE  30 
YEARS   IN  PERIOD  30 

VALUE  OF 

IN] 

[TIAL  GROWING  STOCK- 

-*  60992.01 

VALUES 

DI 

SCOUNTED  TO  PRESENT 

( DOLLARS ) 

COMPOUND 

RATE 
(PERCENT) 

FUTURE 
GROWING 
STOCK 

ALL 
INCOMES 

STOCK 
PLUS 
INCOMES 

ALL 
COSTS 

NET 
PRESENT 
WORTH 

1.0 

38247.39 

52747.65 

90995.04 

15614. 36 

14388. 67 

1.5 

32980.84 

48739.60 

81720.44 

14470. 52 

6257.91 

2.0 

28460. 19 

45142.72 

73602.91 

13440.43 

-829. 53 

2.5 

24576. 89 

41909.46 

66486.35 

12511. 14 

-7016. 80 

3.0 

21238.62 

38998.17 

60236.79 

11671.31 

-12426. 53 

3.5 

18366.76 

36372.40 

54739.15 

10910.99 

-17163.85 

4.0 

15894.35 

34000.11 

49894.46 

10221.45 

-21319. 00 

4.5 

13764.30 

31853.20 

45617.51 

9595.00 

-24969. 50 

5.0 

11927.90 

29906.94 

41834.84 

9024.86 

-28182.04 

5.5 

10343.54 

28139.54 

38483.08 

8505.07 

-31014.00 

6.0 

8975.67 

26531.81 

35507.48 

8030. 34 

-33514.87 

6.5 

7793.90 

25066.81 

32860.71 

7596.04 

-35727. 34 

7.0 

6772.20 

23729.57 

30501.77 

7198.01 

-37688.26 

7.5 

5888.29 

22506.84 

28395.13 

6832.62 

-39429.50 

8.0 

5123.07 

21386.90 

26509.97 

6496.61 

-409178.65 

8.5 

4460. 16 

20359.34 

24819.50 

6187.10 

-42359.61 

9.0 

3885.51 

19414.93 

23300.44 

5901.53 

-43593.10 

9.5 

3387.03 

18545.45 

21932.48 

5637.60 

-44697.13 

10.0 

2954.35 

17743.61 

20697.96 

5393.28 

-45687.33 

10.5 

2578.54 

17002.89 

19581.44 

5166.75 

-46577.33 
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PAGE  TYPE  7  COMPARISON  OF  ALTERNATIVES 

BATCH     TEST  PROBLEM 
TEST  1 

MANAGED, THINNED  AGE  30 
COLUMN  10 


4530 
4621 
4716 
4807 
4897 


4518 
4596 
4677 
4793 
4887 


4992 
5083 
5173 
5273 
5364 


4970 
5068 
5166 
5266 
5357 


5364 
6632 
7893 
0 
0 


5357 
6672 
7831 
0 
0 


90 
100 


110 
120 
130 
140 
150 


COMPARISON  OF  ALTERNATIVES 
BATCH     TEST  PROBLEM 
TEST  1 

MANAGED, THINNED  AGE  30 
COLUMN  40 


65346 
78294 
75197 
61083 
65246 


65182 
77792 
74586 
61283 
65465 


81284 
57210 
54890 
60757 
64798 


81255 
57628 
55403 
61269 
65310 


64798 
80447 
94259 


65310 
81931 
95470 
0 
0 


90 
100 


110 
120 
130 
140 
150 


Apiculture  —  CSU,  Ft.  Colli™ 
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